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1.  INTRODUCTION 


1 

n 


The  earth's  radiation  belts  consist  of  energetic  electrons  and 
ions  (mainly  protons)  which  execute  quasi-periodic  trajectories  under 
the  constraining  influence  of  the  geomagnetic  field.  The  adiabatic  theory 
of  charged-particle  motion  (Northrop  and  Teller,  1960)  provides  the 
kinematical  framework  for  characterizing  radiation-belt  particles.  The 
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kinematical  state  of  a particle  is  defined  by  specifying  the  three  invar-  f 

iant  action  integrals  Jj  and  three  conjugate  phases  associated  with  the 

particle's  gyration  about  a guiding  field  line  (i  = 1),  bounce  motion  | 

between  magnetic  mirror  points  (i  = 2),  and  azimuthal  drift  around  the 

earth  (i  = 3).  This  hierarchy  typically  applies  for  particle  momenta  i 

p«60(Z/L^)  GeV/c,  where  Z is  the  charge  number  (Z  = 1 for  protons 


and  electrons,  Z = 2 for  alpha  particles,  etc.  ).  Moreover,  the  bounce 
and  drift  integrals  are  well-defined  only  for  values  of  the  shell  parameter 
L < 10,  since  the  outer  magnetosphere  is  too  distorted  to  support  adia- 
batic motion. 

Dynamically  interesting  phenomena  are  those  that  violate  the 

adiabatic  invariants,  which  are  respectively  proportional  to  the  action 

2 

integrals  Jj.  The  first  invariant  M = jqj  (Jj/2iTmQc)  = p^^  /2mQB  is 
equal  to  m/m^  times  the  particle's  magnetic  moment.  The  second  invar- 
iant J = J2  = ^P||ds  is  evaluated  along  the  particle's  guiding  field  line, 
and  the  third  invariant  ® = (c/q)Jj  is  equal  to  the  magnetic  flux  enclosed 
by  the  drift  shell.  In  accordance  with  Hamilton- Jacobi  theory,  the  char- 
acteristic frequency  associated  with  the  action  integral  is  given 


Preceding  pege  blank 


-9- 


by  = 3W/9J^,  where  W is  the  particle's  energy  (kinetic  plus 

potential).  Representative  frequency  contours  are  shown  in  Figure  1 
for  particles  mirroring  at  the  equator  of  an  idealirsd  geomagnetic  dipole 
field,  with  no  electrostatic  fields  superimposed  (Schulz  and  Lanzerotti, 
1974). 

Processes  that  violate  either  or  both  of  the  first  two  invariants 
lead  to  particle  diffusion  in  pitch  angle  and/or  energy.  Processes  that 
violate  the  third  invariant  lead  to  radial  diffusion.  In  particular,  processes 
that  conserve  the  first  two  invariants  while  violating  the  third  produce 
a kind  of  radial  diffusion  that  energizes  particles  as  they  diffuse  into  the 
magnetospheric  interior  from  an  external  source,  e.g.  , from  the  plasma 
sheet  or  the  solar  wind  (see  Figure  2). 

The  plasma  sheet  is  an  electrically  resistive  medium  extending 
across  the  magnetospheric  tail  in  a dawn-to-dusk  electric  field  ~1  V/km. 
The  resistivity  is  not  collisional  but  anomalous  (plasma-kinetic),  and  the 
plasma-sheet  temperature  amounts  to  several  keV  (with  T^  > T^).  The 
main  source  of  particles  for  the  plasma  sheet  is  most  likely  the  shocked 
solar  wind,  but  (according  to  Axford,  1970)  plasma  from  the  earth's 
polar  ionosphere  may  also  contribute  via  the  "polar  wind",  which  flows 
out  along  open  field  lines  (see  Figure  3). 

Magnetospheric  plasma  tends  to  drift  as  illustrated  in  Figure  2 
under  the  influence  of  electrostatic  and  magnetostatic  fields.  The  earth's 
rotation  leads  to  a quasi- stable  capture  of  equatorial  and  mid-latitude 
ionospheric  plasma  within  a sharply  bounded  region  known  as  the  plasma - 
sphere.  Here  the  plasma  temperature  is  ~1  eV  or  less,  but  the  equatorial 


-10- 


antours  of  constant  adiabatic  gyration,  bounce,  ajid  drift 
eouency  for  equatorially  mirroring  particles  ingeo- 

agnetic  dipole  field.  The  approximation  of  an  ^diajjatic 

erarchy  of  motions  is  seen  to  fail  around  the  upper 
.rn«rfl  (E  ~ 1 GeV.  L - 8),  where  0.  ~ n,  ~ S' 


PLASMA 

PAUSE 


Fig.  2.  Distribution  (shaded  areas)  and  flow  pattern  (solid  arrows) 
of  magnetospheric  plasma  in  the  equatorial  plane.  Dashed 
arrows  indicate  electric -current  pattern.  The  plasma - 
sphere  is  the  shaded  region  inside  the  plasmapause.  Field 
lines  beyond  the  magnetic  cusp  are  "open",  forming  a 
neutral -sheet  configuration  maintained  by  plasma -sheet 
currents.  However,  the  plasma  sheet  also  extends 
(presumably  by  diffusion)  onto  closed  field  lines  sunward 
of  the  cusp. 


-12- 


MAGNETO- 

SHEATH 


SOLAR 


PLASMA  SHEET 


WIND 


“BOW 

SHOCK 


MAGNETO- 
PAUSE 


F 


MAGNETO- 

SHEATH 


MAGNETOPAUSE 


Fig.  3.  Magn-.io spheric  configuration  in  the  noon-midnight 
meridional  plane.  Shaded  regions  correspond  to 
electron  radiation  belts  (omnidirectional  flux  14^^  2 10^ 
cm"^.sec"*  at  .dnetic  energies  E 2 0.  5 MeV).  Quasi- 
trapping  regions  (bounded  by  dashed  curves)  contain 
the  mirror  points  of  particles  which,  because  of  the 
distorted  magnetic -field  configuration,  cannot 
complete  an  adiabatic  drift  orbit  around  the  earth. 
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density  is  .high  (~  10  e cm"  , according  t''  Chappell  et  al.  , 1970) 

_3 

compared  with  that  in  regions  beyond  the  p'asmapause  (Si  cm  ). 

Since  the  cold-plasma  density  strongly  affects  dispersion  relations 
and  resonance  conditions,  the  plasmapause  (see  Figure  2)  is  also  a 
boundary  for  phenomena  dependent  on  magnetospheric  wave -particle 
interactions. 

In  addition  to  the  above -described  plasma  she^^t,  there  are  other 
important  sources  of  geomagneticaily  trapped  radiation.  Energetic 
solar-flare  protons  (E  ~ 1 — 100  MeV)  can  enter  the  magnetosphere 
more-or-less  directly  from  interplanetary  space.  Such  particles  are 
not  strongly  affected  by  traversal  of  the  magnetosheath,  nor  by  large- 
scale  electric  fields  in  the  magnetosphere.  They  can  enter  the  magneto- 
sphere through  the  tail  (perhaps  also  at  the  day  side  neutral  points)  and 
can  gain  access  to  closed  drift  shells  through  an  intrinsic  breakdown  of 
the  adiabatic  hierarchy,  a breakdown  that  is  characteristic  of  high- 
rigidity  particles  in  the  outer  magnetosphere. 

Solar-flare  protons  striking  the  atmosphere  can  eject  energetic 
neutrons  by  spallation  of  nitrogen  and  oxygen  nuclei.  The  beta  decay 
of  such  albedo  neutrons  within  the  magnetosphere  (a  process  known  as 
SPAND,  for  solar-proton-albedo-neutron  decay)  constitutes  a minor 
source  of  radiation -belt  protons.  The  analogous  reaction,  induced  by 
more  energetic  cosmic-ray  particles  from  the  galaxy,  is  know  as  GRAND 
(for  cosmic-ray-albedo-neutron  decay)  and  accounts  in  a major  way  for 
the  injection  of  energetic  protons  into  the  inner  radiation  zone  (LS  3).  A 
final  source  of  geomagneticaily  trapped  radiation,  operative  during  the 
years  1958  — 62,  consisted  of  the  fission  debris  from  high-altitude  (2  200  km) 
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nuclear  detonations  (e.  g.  . Walt,  1971).  The  beta  decay  of  such  fission 
fragments  led  to  the  formation  of  artificial  radiation  belts  having  various 
widths  in  L and  various  distributions  in  pitch  angle,  and  the  analysis  of 
data  on  such  artificial  radiation  belts  has  provided  valuable  diagnostic 
information  about  the  magnetospheric  environment. 

The  purpose  of  the  present  review  is  to  cover  the  major  develop- 
ments in  radiation-belt  phenomenology  of  the  past  four  years  (1970-1973). 
The  selection  of  topics  is  somewhat  subjective,  and  the  omission  of  some 
possibly  important  contributions  (e.  g. , those  not  understood  by  the  present 
reviewer)  is  unavoidable.  The  intent  is  to  provide  a coherent  and  some- 
what critical  review  of  recent  progress  in  radiation-belt  physics,  and  to 
Suggest  some  promising  ideas  for  further  investigation. 


2.  ADIABATIC  DRIFT  SHELLS 

Although  the  dynamically  interesting  phenomena  of  radiation-belt 
physics  require  (by  the  present  definition)  a violation  of  one  or  more 
adiabatic  invariants,  there  exist  two  interesting  kinematical  phenomena 
(known  as  shell  splitting  and  quasi -trapping)  which  arise  from  the 
magnetosphere's  inherent  azimuthal  asymmetry.  If  the  magnetosphere 
were  azimuthally  symmetric  about  the  earth's  dipole  axis,  the  tr.^.cing  of 
particle  drift  shells  would  consist  simply  of  rotating  each  field  line  about 
the  axis  of  symmetry.  Two  particles  mirroring  on  the  same  field  line. 
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regardless  of  their  energies  or  equatorial  pitch  angles,  would  generate 
the  same  adiabatic  drift  shell. 

However,  the  real  magnetosphere  is  not  azimuthally  symmetric. 
There  is  a day -night  asymmetry  with  respect  to  the  magnetic  field,  and 
a dawn-dusk  asymmetry  with  respect  to  the  magnetospheric  electric  field. 

In  consequence  of  these  asymmetries,  particle  drift  shells  tend  to  inter- 
sect rather  than  coincide.  Two  particles,  even  if  they  have  the  same 
kinetic  energy  on  a given  field  line,  will  (in  general)  execute  different 
drift  shells  unless  they  also  have  the  same  equatorial  pitch  angle  on  the 
j given  field  line.  This  effect  is  known  as  shell  splitting. 

I Evidence  for  shell  splitting  in  a given  field  model  is  obtained  by 

i 

considering  the  drift  trajectory  of  a particle  having  a vanishing  second 
invariant  J.  This  trajectory  must,  of  course,  conserve  the  first  invariant 

I 

i M and  total  (kinetic  plus  potential)  energy  W.  It  must  also  track  a local 

I 
I 

i minimum  in  the  field  magnitude  B (with  respect  to  arc  length  s)  on  each 

field  line  that  it  crosses.  In  analytical  terms,  the  surface  generated 

by  the  drift  paths  of  particles  having  J = 0 is  the  surface  on  which  9B/9s  = 0 
2 2 

9 B/9s  2 0,  where  9/9s  = B • V . It  was  noted  by  Mead  (1964)  that,  in 
his  field  model,  the  J = 0 surface  bifucates  into  two  separate  sheets  on 
the  day  side  of  the  magnetosphere.  Shabansky  and  Antonova  (1968) 

I suggested  that  a particle  having  JssO  will  select  either  the  northern  or 

the  southern  sheet,  according  to  the  instantaneous  value  of  its  bounce 
phase  If  2 bifurcation. 

[■ 
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The  drift  path  of  a particle  having  J = 0 is  a contour  conserving 
both  M and  W on  the  surface  described  above.  (Neglecting  electrostatic 


potentials,  it  is  a path  of  constant  B. ) Shell  splitting  is  indicated  by  an 

2 2 

azimuthal  variation  of  0 B/3s  along  this  path  of  constant  M and  W 
(Roederer  and  Schulz,  1969;  1971). 

It  may  happen  that  a particle  has  well-defined  M and  J on  a given 
field  line,  but  cannot  execute  a complete  adiabatic  drift  shell  without 
encountering  the  magnetospheric  surface  (either  at  the  magnetopause  or 
at  the  neutral  sheet,  Figure  3).  Such  a particle  is  said  to  be  quasi-trapped, 
and  its  mirror  points  are  said  to  occupy  a quasi -trapping  region  for  the 
energy  (W)  in  question.  One  example  of  quasi -trapping  is  the  fate  of 
particles  having  vanishing  kinetic  energy  beyond  the  plasmapause  (see 
Figure  2).  Such  particles  encounter  the  magnetopause  before  completing 
a drift  orbit  of  2ir  radians.  (Despite  the  azimuthal  asymmetry,  there  is 
no  shell  splitting  in  the  limit  of  vanishing  M and  J.  Such  particles 
follow  electrostatic  equipotentials,  which  coincide  with  magnetic  field 
lines  in  most  accepted  models. ) 

The  geometry  of  the  magnetic  field  alone  is  relevant  to  the  motion 
of  those  particles  having  sufficiently  high  energy  to  justify  the  neglect  of 
electrostatic  potentials,  but  having  sufficiently  low  rigidity  to  move 
adiabatic  ally.  The  criterion  on  rigidity  (cf.  Taylor  and  Hastie,  1971) 
can  be  made  more  precise  by  introducing  the  traceless  tensor 


€ 


(pc/qB^)  VB 
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and  requiring  that  each  rectangular  con;ponent  of  c be  mucli  less  than 

unity  in  absolute  value.  This  criterion  is  fortunately  satisfied  by  mosi. 

of  the  particles  one  chooses  to  study  in  the  context  of  radiation-belt 

physics.  Special  methods  of  numerical  analysis,  designed  to  trace  the 

details  of  each  particle  trajectory,  are  obligatory  in  situations  where 

the  components  of  t are  nci  sufficiently  small  (e.  g.  , Smart  et  al.  , 1969). 

2SS 

The  tracing  of  charged-particle  orbits,  whether  by  numerical 
trajectory  analysis  or  by  the  theory  of  adiabatic  invariance,  requires 
the  specification  of  a field  model.  Recent  developmentr  in  field  modeling 
are  reviewed  elsewhere  (e.  g.  , Morfill,  19''5).  It  will  be  sufficient  here 
to  recall  three  rather  old  models.  The  first,  due  to  Mead  (1964),  is  a 
13-term  expansion  of  the  geomagnetic  scalar  potential  V (such  that  B = 

- V V)  in  spherical  harmonics  having  the  north-south  and  dawn-dusk  sym- 

mm 

metries  appropriate  to  the  case  of  a solar  wind  normally  incident  on  a 
dipolar  magnetic  field.  Mead  required  the  solar-wind  plasma  to  undergo 
specular  reflection  from  the  self-consistent  magnetopause  that  forms 
thereby  (Mead  and  Beard,  1964).  The  result  is  a field  pattern  illustrated 


in  Figure  4a. 

If  one  adopts  a microscopic  view  and  allows  the  solar-wind 
plasma  to  penetrate  the  magnetopause  (e.  g,  , Bird  and  Beard,  1972), 
the  resulting  gyration  and  drift  currents  are  found  to  generate  a self- 
consistent  magnetotail  (cf.  Figure  3).  Williams  and  Mead  (1965)  simu- 
lated this  effect  by  adding  to  the  13 -term  model  an  equatorial  current 


-18- 


Fig.  4.  Schematic  representation  of  meridional  field  lines  in  (a) 
the  13 -term  and  (b)  the  3 -term  Mead  magnetospheres 
(solid  curves).  Corresponding  dipole  field  lines  (dashed 
curves)  are  shown  for  invariant  latitudes  A = 65°,  70°, 
75°,  80°,  85°  and  90°,  but  omitted  for  A = 60°.  The 
symbol  X marks  the  location  of  the  nightside  neutral 
line  that  automatically  appears  in  the  3 -term  model 
(Schulz  and  Lanzerotti,  1974). 


strip  across  the  nightside  magnetosphere  from  dawn  to  dusk.  However, 
one  gets  a similar  effect  without  the  current  strip  (see  Figure  4b)  by 
simply  truncating  the  expansion  of  Mead  (1964)  aiter  the  third  term.  The 
resulting  three-term  model,  actually  suggested  by  Mead  (1964),  contains 
the  Sblient  features  of  the  actual  magnetospheric  topology.  It  is  derivable 
from  the  scalar  potential 

I 

I 

V = - (a/r)^z  - g°  (a/b)^z  - sH  (a/b)^  (x/b)  z,  (2) 

where  r is  the  geocentric  distance  (radial  coordinate),  z is  the  distance 
north  of  the  equator  (along  the  dipole  axis),  and  x is  the  projected  distance 
from  the  dipole  axis  (measured  along  the  earth- sun  line,  positive  on  the 
night  side).  The  expansion  coefficients  and  g™  are  Schmidt- 

normalized  (Chapman  and  Bartels,  1940).  The  symbols  ^ and  b 
represent  the  radius  of  the  earth  and  the  geocentric  "stand-off"  distance 
to  the  subsolar  point  on  the  magnetopause,  respectively. 

Using  the  field  model  of  Williams  and  Mead  (1965),  Roederer  (1967) 
has  studied  the  qua  si -trapping  phenomenon  in  situations  for  which  elec- 
trostatic potentials  can  legitimately  be  neglected.  He  has  found  that  a 
particle  having  one  mirror  point  in  each  of  the  two  qua:;: -trapping  regions 
on  the  noon  meridian  (regions  Q in  Figure  3)  will  fail  to  complete  a drift 
orbit  of  2tt  radians,  conserving  M and  J.  Such  a particle  would  have  to 
preserve  the  magnitude  of  its  mirror  field  B in  the  course  of  azimuthal 
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2 2 

drift,  while  conserving  also  the  quantity  K (s  J /SrriQM).  This  proves 

to  be  impossible,  and  the  particle's  drift  shell  instead  crosses  the 

nightside  boundary  between  closed  and  open  field  lines. 

The  foregoing  conclusion  applies  to  bounce  trajectories  that  cross 

the  equatorial  plane,  i.e.  , to  mirror-field  intensities  B that  exceed  the 

" " m 

magnitude  of  B in  the  equatorial  plane  along  a given  field  line,  in  the  noon 
meridian.  Roederer  (1967)  did  not  consider  the  motion  of  particles 
trapped  in  the  higher -latitude  field  minima  that  exist  near  the  magneto- 
pause on  the  day  side.  Thus,  his  discovery  of  a nightside  qviasi-trapping 
region  (marked  Q in  Figure  3),  characteristic  of  particles  mirroring  too 
near  the  magnetic  equator,  must  be  received  with  caution  (cf.  Shabansky 
and  Antonova,  1968). 

Particles  having  K = 0 do,  however,  show  evidence  of  quasi- 
crapping in  tht  three-term  field  model  of  Mead  (1964).  It  can  be  shown 

2 2 

(cf.  Roederer  and  Schulz,  1969)  that  the  value  of  9 B/9s  is  correctly 
given,  in  the  equatorial  plane  of  this  field  model,  by  the  expression 


r^B(a^B/9s^) 


9B^  + 9g°(a/b)^[3B  + 2g°(a/b)^] 

+ 6(g2S^(a/b)^(r/b)^(  1 + 15  cos^^P) 

- *v/3  g^  (a/b)^(r /b)  [39B  + 48  g|^(a/b)^]cos  (p. 


where  if  is  the  longitude  measured  eastward  from  t'- e midnight  meridian. 

2 2 

The  azimuthal  variation  of  a positive  (9  B/9s  )^^  demonstrates  the  exis- 
tence of  shell  splitting  in  this  model.  Moreover,  it  follows  from  (3) 

2 2 

thet  9 B/9s  vanishes  in  the  equatorial  plane  along  the  contour  given  by 
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^J3  cos  if  = (2g2^7g°)(r/b)'^  + (3g  °/7g2  )(b/r)^ 

- {3gf/7g2)(b/r).  (4) 

This  contour  marks  the  aforementioned  bifurcation  of  the  K = 0 surface. 
Moreover,  the  magnitude  of  B vanishes  in  the  equatorial  plane  along  the 
contour  given  by 

s!i  cos  if  = - ig^/g^)  (b/r)^  - (gj^/g^)  (b/r).  (5) 

This  contour  marks  the  boundary  between  closed  and  open  field  lines 

(see  Figure  4b).  Both  of  these  important  contours  are  plotted  in  Figure  5, 

using 


g®  = 1.  24  g®  = - 2.  56  g2  = - 0.  31  gauss,  (6) 

as  recommended  by  Mead  (1964).  Both  contours  intersect  the  equatorial 

trace  of  the  magnetopause  calculated  self-consistently  (i.  e.  , in  the  model 

containing  many  more  than  three  spherical-harmonic  coefficients)  by 

Mead  and  Beard  (1964).  The  analytical  continuation  of  B = 0 and 
2 2 

(a  B/9s  )^^  = 0 beyond  the  Mead-Beard  surface  is  indicated  by  dashed 
curves,  but  presumably  has  no  physical  meaning. 

It  is  convenient,  following  Stone  (1963),  to  label  a field  line  with 
the  parameter 


(7) 
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2 2 

Fig,  5.  Equatorial -plane  contours  of  constant  B and  of  9 B/9s  =0 
in  the  three -term  Mead  field.  Dashed  curves  represent 
the  analytical  continuation  of  contours  B = 0 and  9^B/9s^  = 0 


beyond  the  magnetospheric  surface  calculated  by  Mead  and 
Beard  (1964).  Dotted  curves  represent  azimuthal -drift 
trajectories  of  particles  having  J = 0.  The  contour 
Lq  = 1.2b/a  illustrates  quasi -trapping  (cf.  Figure  3),  and 
the  contour  Lq  = 0.  8b/a  illustrates  dayside  bifurcation  of 
the  azimuthal -drift  trajectory. 


where  is  the  minimum  value  of  B along  the  field  line.  Thv 

parameter  Lq  thus  serves  also  as  a label  for  the  drift  path  of  a particle 

having  K = 0.  Figure  5 thus  illustrates  (by  means  of  dotted  curves)  both 

stable  trapping  (Lq  = 0.  4 b/a)  and  quasi-trapping  (Lq  = 1 . 2 b/a)  at  K = 0, 

as  well  as  an  intermediate  case  (Lq  =0.8  b/a)  in  which  the  drift  path 
2 2 

bifurcates  at  9 B/9s  = 0.  Further  analysis  would  be  required  to  determine 
whether  the  drift  path  Lq  = 0.  8 b/a  closes  within  the  magnetosphere  in  this 
model  of  the  earth's  field. 

The  above  results  tend  to  confirm  the  suggestion  by  Roederer 
(1967)  that  there  is  a nightside  quasi -trapping  region  containing  the 
mirror  points  of  particles  having  small  values  of  K.  In  both  studies, 
however,  the  field  models  used  contain  elements  of  inconsistency.  In 
the  present  study,  the  three-term  representation  of  the  field  is  manifestly 
inconsistent  with  the  magnetopause  determination  of  Mead  and  Beard  (1964), 
which  would  require  at  least  the  IT*-  ' • ' field  model  for  self-consistency, 

in  the  case  of  Roederer 's  study,  thv  • u model  contained  an  infinitely 
wide  current  sheet  (Williams  a.nd  Mead,  1965)  extending  from  the  dawn 
horizon  to  the  dusk  horizon  The  currents  in  that  model  thus  did  not 
properly  close  over  the  rounded  surface  of  the  magnetotail  (cf.  Figure  3). 

It  remains  an  open  question  whether  any  particles 

having  K = 0 can  show  quasi-trapping  in  a self-consistent  model  of  the 
geomagnetic  field  and  its  boundary.  A study  based  on  the  13 -term  model 
of  Mead  (1964)  would  help  to  resolve  this  essentially  topological  question, 
although  the  inclusion  of  a self-consistently  modeled  neutral  sheet  would 


be  necessary  in  order  to  resolve  the  question  of  quasi -trapping  in  the 
real  geomagnetic  field. 

Although  presented  above  as  a means  of  particle  escape  from  the 
magnetosphere,  quasi-trapping  can  equally  well  be  viewed  as  a means  of 
particle  access.  The  equations  of  motion  a -e  quite  reversible  in  space 
and  time.  However,  such  access  from  outside  the  magnetosphere  would 
be  only  transient  if  unaccompanied  by  some  means  of  radial  diffusion 
from  the  partial  (quasi-trapped)  drift  shell  to  a closed  drift  shell.  One 
may  consider  in  this  context  the  access  of  solar-wind  plasma  to  the 
earth's  plasma  sheet,  by  way  of  the  dayside  polar  "cusp"  (e.  g.  , Frank, 
1971),  which  extends  roughly  from  either  neutral  point  down  to  the  earth's 
surface  (see  Figure  3). 

It  has  been  noted  above,  following  equation  (3),  that  shell  splitting 
is  present  even  under  the  simpler*  idealization  of  the  day-night  magneto- 
spheric  asymmetry.  Shell  splitting  is,  of  course,  present  also  in  more 
nearly  realistic  field  representations,  such  as  the  model  of  Williams  and 
Mead  (1965).  Roederer  (1967)  traced  the  adiabatic  drift  shells  of  this 
field  model  and  thereby  estimated  the  amount  by  which  drift  shells  are 
split  apart.  Pfitzer  et  al.  (1969)  carried  the  analysis  a significant  step 
further  by  using  a shell-tracing  program  to  map  the  distribution  function 
of  geomagnetically  trapped  particles  from  one  longitude  to  another.  They 
used  data  from  an  elliptically  orbiting  spacecraft  (CXlO-3)  to  determine 
the  radial  distribution  of  electron  intensity  at  each  equatorial  pitch  angle. 
Applying  Liouville's  theorem,  they  then  used  the  shell-tracing  program 
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to  predict  the  diurnal  variations  of  electron  intensity  at  synchronous 
altitude  (specifically  at  the  site  of  ATS  1).  The  good  agreement  between 
prediction  and  reality  is  illustrated  in  Figure  6,  where  x denotes  the 
cosine  of  the  equatorial  pitch  angle,  is  the  tail-field  intensity  (an 

_5 

adjustable  parameter  of  the  Mead-Williams  model),  and  1 Y = 10  gauss. 

It  is  interesting  to  notice  that  shell  splitting  thus  leads  to  an 
anomalous  pitch-angle  distribution,  i.  e.  , a distribution  having  a relative 
minimum  at  x = 0,  in  the  general  vicinity  of  midnight.  The  pitch -angle 
distribution  is  normal  (essentially  compatible  with  the  lowest  eigenmode 
of  pitch-angle  diffusion)  at  the  noon  meridian  (1200  i^T).  The  explanation 
is  not  difficult.  At  L2  5,  che  distribution  of  electron  intensity  at  fixed 
energy  decreases  with  increasing  drift-shell  radius.  In  other  words,  the 
radiation  resides  in  a belt.  Particles  having  x = 0 travel  a path  of 
constant  equatorial  B,  and  must  pass  nearer  the  earth  at  midnight  than 
at  noon  in  order  to  do  so.  Particles  having  x = 0.  42  (for  example) 
mirror  in  the  relatively  less  distorted  inner  magnetosphere,  and  so 
execute  a more  nearly  symmetrical  drift  shell.  Thus,  the  ATS-1  space- 
craft at  midnight  samples  electrons  from  a larger  drift  shell  for  x = 0 
than  for  x = 0.  42;  the  ordering  of  shell  diameters  is  reversed  when  the 
spacecraft  is  at  noon. 

One  immediately  observes  from  Figure  6 that  the  midnight  pitch- 
angle  distribution  is  more  strongly  anomalous  in  the  higher-energy  channel. 
This  empirical  fact  was  scrutinized  more  thoroughly  by  Lucas  and  Brice 
(1973),  using  the  electron  data  of  West  et  al.  (1973).  Representative 
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Fig,  6.  Diurnal  variation  of  electron  fluxes  and  magnetic  field 
(data  points)  observed  at  synchroncas  altitude  (ATS  1), 
as  compared  with  predictions  (dashed  and  solid  curves) 
based  on  Mead -Williams  field  models  (b/a,B^)  and  dis- 
tribution function  f(M,  J,®)  deduced  from  OGO -3  electron 
data  (Pfitzer  et  al.  , 1969). 
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nightside  (left  panel)  and  morningside  (right  panel)  pitch-angle  distributions 
are  shown  in  Figure  7.  Even  at  the  lower  energies,  where  the  nightside 
distribution  remains  normal,  one  detects  a relative  flattening  of  the  pitch- 
angle  distribution  on  the  night  side.  West  et  al.  (1973)  obtained  pitch-angle 
distributions  at  a great  variety  of  longitudes.  On  examining  their  charac- 
terizations (here  called  normal  vs  anomalous),  one  detects  a gradual 
rotation  of  the  predominant  asymmetry  from  day-night  toward  dawn-dusk 
with  decreasing  electron  energy.  This  rotation  is  qualitatively  compatible 
with  the  incipient  contribution  of  the  magnetospheric  electric  field  to 
drift -shell  splitting  (see  Roedevt  r and  Schulz,  1971;  Stern,  1971). 

Although  shell  splitting  itself  is  a purely  adiabatic  phenomenon 
(conserving  all  three  invariants),  it  can  lead,  in  the  presence  of  particle 
collisions  with  the  atmosphere,  to  an  enhanced  form  of  radial  diffusion 
(called  "neoclassical"  diffusion  in  the  language  of  laboratory  plasma  physics). 
This  may  be  an  important  piocess  for  inner-zone  electrons  (L^l.2), 
where  other  known  radial-diffusion  mechanisms  are  relatively  ineffective. 

The  shell  splitting  at  such  low  E values  arises  not  much  from  the  gross 
magnetospheric  asymmetries  discussed  above,  but  rather  from  higher- 
order  internal  multipoles  (magnetic  anomalies)  of  the  earth's  main  field. 

The  magnitude  of  this  internal  shell  splitting  is  derivable  from  Figure  8 
(Roederer  et  al.  , 1973),  in  which  the  subscript  0 is  used  to  emphasize 
that  (J  B/'ds  has  been  evaluated  along  a path  of  constant  B ( = Bq)  on 
a warped  surface  defined  by  the  condition  8B/9s  = 0 (see  above).  The 
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Fig.  7,  Typical  nighttime  (18  September  1968)  and  morningside 
(30  March  1968)  pitch-angle  distributions  of  energetic 
electrons  at  L ~ 7 (from  West  et  al.  , 1973).  The  data 
were  transformed  to  eqxiatorial  latitude;  the  dashed 
portions  of  the  curves  are  extrapolations  (Lucas  and 
Brice,  1973). 


A 


minimum-B  surface  is  no  longer  planar,  as  it  was  in  equations  (2)-(5). 
The  definition  of  Lq  is  still  given  by  (7)  with  the  interpretation  that  g^^ 
is  to  be  evaluated  in  offset-dipole  coordinates  (Chapman  and  Bartels, 

1940).  Shell  splitting,  here  as  elsewhere,  is  indicated  by  the  azimuthal 

2 2 

variation  of  d B/9s  along  a path  of  constant  B ( = Bq)  on  the  surface 
K = 0.  The  main  features  discernible  at  L^-v-l  are  the  South  American 
(f^»0°-30°)  and  South  African  {(p  a 40°-110°)  anomalies.  The  main 
asymptotic  contribution  to  shell  splitting  (actually  dwarfed  by  gross  mag- 
netospheric  asymmetries  at  Lq2  2)  arises  from  the  internal  octupole 
(n  = 3,  m = 2). 

3.  PARASITIC  PITCH -ANGLE  DIFFUSION 

Following  a major  disturbance  (either  natural  or  artificial)  of 
the  magnetospheric  environment,  one  generally  observes  a temporal 
decay  of  enhanced  electron  intensities  toward  the  levels  characteristic 
of  magnetically  quiet  periods.  For  kinetic  energies  E Z 100  keV,  the 
decay  is  typically  exponential,  yielding  a lifetime  t that  depends  upon 
E and  L.  A compilation  of  such  electron  lifetimes,  based  on  observa- 
tions of  the  integral  flux  at  E >0.  5 MeV,  is  illustrated  in  Figure  9.  The 
lifetimes  ~ 200  days  (filled  circles)  at  the  uvo  lowest  L values  are 
compatible  with  the  theory  of  atmospheric  scattering  (Walt,  1966). 
However,  the  decay  rates  observed  beyond  L %1.  5 clearly  demand 
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another  explanation.  It  is  by  now  generally  accepted  that  wave -pa rticle 
interactions  are  the  responsible  scattering  agent,  and  that  w thin  the 
plasmasphere  the  relevant  waves  are  electromagnetic. 

Since  the  decay  of  electron  intensities  at  E ^100  keV  is  typically 
exponential,  it  must  be  that  the  wa^'es  responsible  for  the  particle  scat- 
tering are  not  primarily  generated  by  the  particles  being  scattered.  The 
waves  may  have  arisen  from  some  other  source,  e.  g.  , a population  of 
lower-energy^  particles  residing  elsewhere  in  the  magnetosphere,  or 
they  may  have  arisen  in  part  from  earth-based  radio  sources.  A minor 
fraction  of  the  magnetosphe ric  VLF  wave  energy  can  be  attributed  to 
"whistlers"  created  by  lightning  discharges.  Since  the  particles  of 
interest  here  thus  "feed  upon"  wave  energy  not  of  their  own  creation, 
their  pitch-angle  diffusion  is  said  to  be  parasitic.  No  evil  connotation 
is  intended. 

The  empirical  consequence  of  magnetospheric  wave-particle 
interactions  is  a particle  decay  rate  that  increases  monotonically  with 
L.  This  'onclusion  is  reinforced  by  Figure  10,  which  is  a compilation 
of  post-storm  electron  lifetimes  in  three  energy  channels  (Williams  et  al 
1968).  Here  the  lifetime  t is  seen  to  be  a monotonically  increasing 
function  of  particle  energy  and  (as  noted)  a monotonically  decreasing 
function  of  L. 

Early  efforts  to  model  the  interaction  between  gcomagnetically 
trapped  particles  and  electromagnetic  waves  generally  focused  on  an 


equatorial  interaction  via  cyclotron  resonance,  as  specified  by  the 
condition 

w - kjVji  = n,  (8) 

where  w/2tt  is  the  frequency  of  the  wave,  fl/Zir  is  the  relativistic 
gyrofrequency  of  the  particle,  and  k||V|j  is  the  Doppler  shift  associated 
with  the  particle's  component  of  motion  along  B.  The  main  difficulty 
with  (8)  is  that  particles  mirroring  too  near  the  equator  fail  to  interact 
effectively  with  a reasonable  wave  spectrum. 

The  problem  of  parasitic  pitch-angle  diffusion  has  recently 
been  examined  in  greater  generality  by  Lyons  et  al.  71,  1972)  and 
by  Lyons  and  Thorne  (1972).  They  considered  not  only  the  case  of 
wave  propagation  (j^  parallel  to  as  covered  by  (8),  but  more  gen- 
erally the  case  of  wave  propagation  oblique  to  This  expanded 
outlook  leads  tr>  resonance  under  conditions  for  which 

w - k||V||  = nn,  (9) 

wlu:rc  n is  any  integer.  The  case  n = 1 recovers  the  primary  cyclo- 
tron resonance  in  the  form  of  (8),  but  additional  cyclotron  resonances 
(n  = 2,  3,  4, ... ) facilitate  the  interaction  of  higher-energy  electrons  with 
the  wave  spectrum.  Moreover,  the  Landau  r“sonance  (n  = 0)  enables 
particles  mirroring  very  near  the  equator  to  interact  effectively  with  a 
reasonable  wave  spectrum  (see  Figure  11). 


DIFFUSION  COEFFICIENT  (arbitrary  units) 
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L=4,  BOUNCE  AVERAGED 


30"  60" 

EQUATORIAL  PITCH  ANGLE 


Bounce -orbit  averaged  contributions  from  cyclotron 
resonance  (including  harmonics)  and  Landau  resonance 
to  the  pitch -angle  diffusion  coefficient  for  various 
electron  energies  at  L = 4,  as  determined  by  Lyons 
et  al.  (1972). 
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Recognising  ttu^t  wave -particle  interactions  need  not  be  confined 
to  the  equatorial  plane,  Lyons  et  al.  (1972)  assumed  a uniform  distri- 
bution of  wave  energy  throughout  the  plasmasphere  and  performed  a 
properly  weighted  bounce  average  to  obtain  • the  diffusion 

coefficient  for  the  cosine  of  the  equatorial  pitch  angle.  They  considered 
a broad  Gaussian  wave  spectrum  (half-width  = 300  Hz  at  1/e)  centered  at 
600  Hz,  and  imposed  a lower  cut-off  at  300  H^.  The  angular  distribution 
of  wave  vectors  was  taken  as  proportional  to  the  function 


0(k)  = exp  [0.04  - 0.04  (k.  B)‘^]. 


(10) 


Lyons  et  al.  (1972)  tentatively  neglected  the  electric-field  component 
of  the  waves,  and  so  assured  that  the  resulting  pitch-angle  diffusion  would 
be  perfectly  elastic.  This  assumption  justified  a diffusion  equation  (for 
the  phase-space  distribution  function  f)  of  the  form 


8£  . 1 d 

TF  ■ xT(y)  ”5x 


[■‘T(y>D„  -II-] 


(11) 


E,L 


which  Lyons  et  al.  (1972)  solved  for  its  longest-lived  eigenfunction  gQ(x). 
The  Jacobian  xT(y)  is  well  approximated  by  uce  of  the  formula 


1/2, 


T(y)  = 1.3802  - 0.3198  (y+ y') 
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(12) 


2 1/2 

given  by  Lenchek  et  al.  (1961),  where  y h (1  - x ) . The  eigenfunc- 

tions gjj(x)  are  required  to  be  even  functions  of  x and  to  vanish  at  the 
edge  of  the  loss  cone  (shaded  area  in  Figure  11). 

Inserting  a root-mean- square  wave  field  of  35  mY,  Lyons  et  al. 
(1972)  obtained  the  lifeames  shown  in  Figure  12.  The  corresponding 
eigenfunctions  gg(x)  are  shown  in  Figure  13,  superimposed  with  arbi- 
trary normalization  on  a collection  of  data  provided  by  H.  I.  West 
(personal  communication).  The  "shoulders"  of  the  predicted  pitch-angle 
distribution  in  Figure  13  correspond  to  the  "bottleneck"  in  the  pitch-angle 
diffusion  coefficient  in  Figure  11.  It  is  a remarkable  triumph  of  such  a 
simply  formulated  theory  to  find  that  these  "shoulders"  in  the  predicted 
distribution  function  correspond  so  closely  in  pitch  angle  to  similar 
j features  in  the  observational  data.  (R.  W.  Ficdricks  points  out,  in  a 

i personal  communication,  that  the  integration  of  Lyont'  pitch-angle 

distribution  over  the  angular  aperture  of  West's  instrument  would  further 
improve  the  quantitative  agrtument.  ) 

Inclusion  of  the  wavelike  electric  field  in  such  an  analysis  is 
straightforward  in  principle.  The  wavo -particle  interaction  in  this  case 
I becomes  slightly  inelastic  at  high  energies  (E  d 100  keV),  and  rather 

I strongly  inelastic  at  low  energies  (E  -50  keV).  Formally  speaking,  the 

I diffusion  coefficient  becomes  a second- rank  tensor  accommodating 

I stochastic  changes  in  E and  x,  as  well  as  their  correlation.  The 

procedure  for  evaluating  this  diffusion  tensor  has  been  outlincu  by  Lyons 
I (1974). 
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region  during  the  decay  phase  toilowing  an  injection  event.  % 
made  by  Lyons  et  al.  (1972),  observations  by  H,  I,  West  (pe 
tion,  1971). 
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4.  FORMATION  OF  INNER  ELECTRON  BELT 

It  is  underst«adably  customary  to  detect  radiation-belt  particles 
by  means  of  instruments  calibrated  with  respect  to  energy  and  angle  of 
incidence,  rather  than  with  respect  to  the  adiabatic  invariants.  Thus, 
although  the  adiabatic  invariants  would  have  been  the  more  natural 
coordinates  for  describing  the  distribution  of  geomagnetically  trapped 
radiation,  most  available  maps  of  radiation  intensity  provide  flux  pro- 
files plotted  against  L at  fixed  energy  and  equatorial  pitch  angle.  Electron- 
flux  profiles  obtained  in  this  format,  especially  during  geomagnetically 
quiet  periods,  define  two  seemingly  distinct  radiation  belts  separated  by 
a deep  "slot".  This  pattern  is  illustrated  in  the  left  panel  of  Figure  14. 

The  center  of  the  "slot"  is  located  typically  between  L = 2.  5 and  L =3.5, 
depending  on  the  particle  energy  considered. 

The  quiet-time  "slot"  region  becomes  filled  in  with  particles 
during  a magnetic  storm,  presumably  because  of  enhanced  radial  diffusion. 

The  result  is  illustrated  in  the  right  panel  of  Figure  14.  At  the  lower  energies 
(E'^SOO  keV  in  this  example)  the  replenishment  is  so  complete  that  the 
two-zone  structure  qualitatively  disappears.  At  the  higher  energies  the 
replenishment  serves  to  shift  the  apparent  "slot"  position  to  L values 
less  than  L = 2. 

During  the  period  of  decay  following  a storm,  the  electron-flux 
profiles  return  to  their  respective  quiet-time  configurations,  as  repre- 
sented by  the  left  panel  of  Figure  14,  and  this  pattern  remains  essentially 
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stationary  until  the  next  major  geomagnetic  disturbance.  One  could 
imagine  that  the  process  (presumably  pitch ^angle  diffusion)  responsible 
for  the  decay  of  radiation  intensity  is  extinguished  upon  restoration  of 
the  quiet-time  configuration,  and  that  particle  motion  thereafter  remains 
perfectly  adiabatic.  This  would  be  wrong.  A more  reasonable  interpre- 
tation of  radiation-belt  dynamics  is  that  the  processes  responsible  for 
radial  diffusion  and  pitch-angle  diffusion  during  disturbed  periods  persist 
(perhaps  with  modified  intensity)  during  the  intervening  quiet  periods, 
and  that  the  static  flux  profiles  observed  during  quiet  periods  represent 
a detailed  balance  between  the  effects  of  radial  diffusion  and  pitch-angle 
diffusion. 

Lyons  and  Thorne  (1973)  have  explored  this  likely  possibility  by 
solving  the  radial-diffusion  equation 


8.  8L  Ll[r  I-*-  8lJ„_k  "c 


(13) 


for  the  quasi-static  phase-space  distribution  function  f.  which  is  eqxial 
to  J^/2mgMB,  where  denotes  the  differential  electron  flux  per  unit 
solid  angle  in  a direction  perpendicular  to  ^ Lyons  and  Thorne  (1973) 
included  loss  terms  both  for  atmospheric  Coulomb  scattering  (t^)  and 
for  wave -particle  interactions  (t^).  They  derived  the  former  lifetimes 
from  Walt  (1966),  and  the  latter  from  Figure  12.  They  considered  the 
radial  diffusion  to  be  caused  by  electrostatic  impulses  (thereby  defining 


-43- 


an  M-dependent  form  for  and  solved  (13)  for  f at  K = 0.  The 

results,  normalized  to  a common  value  at  L = 5.  5,  are  shown  in  the 
left  panel  of  Figure  15.  Since  f must  vanish  at  L«l,  it  is  inevitable 
in  this  formulation  of  the  problem  that  f will  emerge  as  a monotonic 
function  of  L.  at  fixed  M.  Thus,  the  left  panel  of  Figure  15  contains 
no  hint  of  a two-zone  structure. 

Lyons  and  Thorne  (1973)  next  applied  these  results  for  f at 
fixed  M to  obtain  profiles  of  at  fixed  kinetic  energy  E,  using  an 
observed  energy  spectrum  at  L = 5,  5 to  normalize  the  distribution 
function.  The  outcome  of  this  transformation  is  shown  in  the  right  panel 
of  Figure  15.  The  two-zone  structure  thus  arises  quite  naturally  from 
a monotonic  distribution  in  phase  space,  a distribution  consistent  with 
the  radial  diffusior  of  energetic  electrons  from  an  external  source 
(beyond  L = 5.  5)  to  ,in  atmospheric  sink.  Moreover,  the  "slot"  between 
the  radiation  belts  becomes  increasingly  pronounced  with  increasing 
energy,  in  agreement  with  the  observational  data. 

The  next  reasonable  step  in  this  area  of  investigation  would  be 

to  combine  radial  diffusion  and  pitch-angle  diffusion  explicitly  in  the 

same  equation.  T’'e  effective  lifetimes  t and  t used  in  (13)  cor- 

c w ' 

respond  to  a pitch-angle  distribution  in  its  lowest  eigenmode.  However, 
as  Walt  (1970)  has  pointed  out,  radial  diffusion  can  have  the  effect  of 
distorting  the  pitch-angle  distribution.  When  both  processes  act 
simultaneously,  the  radial  profile  and  pitch-angle  distribution  should 
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Fig.  15.  Predicted  steady-state  profiles  of  J^/MB  at  constant  M (left  panel) 
normalized  to  a common  value  at  L = 5.  5,  and  of  at  constant  E 
(right  panel)  normalized  by  prescribing  fue  energy  spectrvim  at 
L - 5.  5 (from  Lyons  ind  Thorne,  1973). 


really  be  determined  in  a self-consistent  manner,  by  including  both 
radial  diffusion  and  pitch-angle  diffusion  on  an  equal  footing  in  the  basic 
equations. 

5.  ARTIFICIAL  RADIATION  BELTS 

Prior  to  the  international  treaty  of  1963,  both  the  USA  and  the 
USSR  had  conducted  nuclear  tests  in  the  atmosphere  and  in  space.  The 
high-altitude  detonations  (see  Table  1)  were  the  source  of  large  numbers 
cl  energetic  electrons  that  became  trapped  in  the  geomagnetic  field, 
thus  forming  artificial  radiation  belts.  Some  of  the  spacecraft  aloft  at 
the  time  of  the  high-altitude  detonations  suffered  severe  radiation  damage. 
Others  launched  immediately  following  several  of  the  respective  detona- 
tions, however,  served  to  provide  a wealth  of  data  on  the  dynamics  of 
artificial  radiation  belts.  Moreover,  since  the  dynamical  processes  ^ 
affecting  artificial  radiation  belts  after  their  formation  are  essentially 
the  same  processes  that  affect  natural  radiation  belts,  one  thus  acquired 
large  amounts  of  data  concerning  the  behavior  of  radiation  belts  in  general. 
Interpretation  of  these  data  continues,  even  at  the  present  time. 

Several  of  the  artificial  detonations  provided  the  observer  with 
particularly  interesting  distributions  of  geomagnetically  trapped  radia- 
tion. For  example,  the  Soviet  detonation  of  28  October  1962  yielded  an 
anomalous  pitch-angle  distribution,  i.  e.  , an  of.^-equatorial  maximum 
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in  the  radiation  intensity  at  L»1.9  (Roberts,  1969).  The  detonation  of 
1 November  1962  yielded  a radiation  belt  confined  to  an  especially  narrow 
interval  in  L (see  Figure  16).  These  two  space  exper’.msnts  thus  led 
naturally  to  analytical  studies  of  pitch-angle  diffusion  and  radial  diffusion, 
respectively. 

Roberts  (1969)  fitted  the  anomalous  pitch-angle  distribution  following 
the  28  October  event  to  a superposition  of  the  three  lowest  eigenfunctions 
of  (11),  making  the  convenient  approximation  that  T(y)nT(l)  outside  the 
loss  cone.  He  thus  accounted  for  the  disappearance  of  the  off -equatorial 
maximum  in  radiation  intensity  on  a shorter  time  scale  than  the  decay  of 
this  artificial  belt  as  a whole. 

Brown  (1966)  fitted  the  narrow  radial  profile  following  the 


1 November  event  to  a temporally  expanding  Gaussian.  He  found  the  squared 

.5 

full  width  at  half  maximum  (FWHM)  to  increase  by  about  7X10  units  of 
L^per  day.  Schulz  and  Lanzerotti  (1974)  translated  this  result  into  an 
equivalent  radial-diffusion  coefficient  6X  10"^  day”^  at  fixed 

energy.  Newkirk  and  Walt  (1968a)  had  earlier  performed  a more  sophis- 
ticated analysis,  transforming  the  observed  profile  into  an  equivalent 
distribution  function  at  fixed  M and  K (=  0)  by  assuming  an  equilibrium 
fission  spectrum  (see  below)  for  the  initial  profile.  They  similarly 
obtained  6X10  ^ day  but  their  method  of  analysis  was  superior 

since  it  properly  conserved  M and  K under  radial  diffusion.  In  neither 
analysis  could  radial  diffusion  account  for  more  than  about  one -third  of 


RELATIVE  INTENSITIES 


Fig.  16.  Evolution  of  inner -zone  electron -flux  profile 
(E  > 1.9  MeV,  omnidirectional)  observed  on 
Explorer  15  following  high -altitude  nuclear 
explosion  of  1 November  1962  (Brown,  1966). 
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the  observed  decline  in  maximum  electron  intensity  (the  peak  at  L = 1.  765). 

For  thi.3,  one  must  postulate  pitch-angle  diffusion  with  a characteristic 

lifetime  of  15-20  days,  as  in  (13). 

The  mechanism  for  creating  an  artificial  radiation  belt  has  been 

described  by  Davidson  and  Hendrick  (1971).  The  high-altitude  nuclear 

detonation  first  creates  a partially  ionized  plasma  that  expands  hydro- 

dynamically  until  the  plasma  pressure  (P)  only  balances  the  magnetic 
2 

pressure  (B  /8tt).  Plasma  expansion  along  ^ is  unimpeded,  except  per- 
haps by  gravity  (since  the  plasma  ions  are  heavy,  being  fission  products). 
Plasma  expansion  across  B is  presumably  mediated  by  turbulent  dif- 
fusion  until  p(&  8irP/B  ) decreases  to  about  unity.  The  plasma  ions  and 
electrons  thereafter  expand  along  B,  and  the  neutral  fission  products 
expand  irrespective  of  B.  The  fission-product  nuclei  are  undergoing  beta 

MM* 

decay  all  the  while,  and  some  of  the  resulting  (beta)  electrons  are  emitted 
with  velocity  vectors  compatible  with  trapping  by  the  geomagnetic  field. 
Although  injected  over  a relatively  narrow  range  of  longitudes,  these  are 
the  electrons  that  ultimately  constitute  the  artificial  radiation  belt.  The 
belt  forms  by  drift-phase  mixing  (dispersion  in  since  the  constituent 

eiectrons  differ  somewhat  from  each  other  in  energy,  equatorial  pitch 
angle,  and  L value  (and  therefore  in  drift  frequency  n^/Zir).  The  artifi- 
cial radiation  belt  thereafter  becomes  subject  to  the  natural  processes 
leading  to  radial  diffusion  and  pitch-angle  diffusion,  and  so  the  deter- 
ministic history  of  each  individual  particle  becomes  obscured. 
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Tha  solid  curves  in  Figure  17  (based  on  Crowther  and  Harless, 

1971)  represent  the  prediction  of  such  a model  for  a nuclear  detonation, 

2? 

having  a fission  yield  of  one  megaton  ( = 4.186  X 10  " erg),  at  altitude  200  km 
on  the  field  line  L = L . The  ordinate  here  is  a damage -equivalent  flux 
of  1-MeV  electrons,  but  can  be  considered  for  present  purposes  as  an 
integral  flux  of  the  fission  electrons.  The  form  of  the  spectrum,  following 
Carter  et  al.  (1959),  should  be  well  approximated  by  the  function 


J^^(E)  « (v/c)  exp  [-  0.  2938  (Y  - 1)  - 0.  0144  (Y  - i)^],  (14) 


where  v is  the  particle  speed,  c is  the  speed  of  light,  and  Y is  the 

ratio  of  relativistic  mass  (m)  to  rest  mass  (m^).  As  a rough  scaling  law 

to  fission  yields  Y other  than  ore  megaton,  the  belt  width  (in  L)  should 

1 /3 

be  proportional  to  Y , and  the  peak  radiation  intensity  should  be  pro- 
portional to  Y^^^. 

10-4 

The  dashed  line  in  Figure  17  represents  a flux  I^^~7  X 10  L 
Z 1 

cm  sec  . This  is  the  maximum  stable  intensity  obtained  from  the 

theory  of  Kennel  and  Petschek  (1966)  for  a pitch-angle  distribution  having 

normal  anisotropy.  A larger  flux  of  electrons  would  lead  to 

the  spontaneous  generation  of  electromagnetic  noise  in  the  whistler  mode, 

and  the  resulting  pitch-angle  diffusion  would  quickly  reduce  the  particle 

intensity  to  • The  stable -trapping  limit  of  Kennel  and  Petschek  (196b) 

stands  far  below  the  condition  P = 1 (solid  line),  under  which  the  particle 

2 

energy  density  equals  B /8ir.  Thus,  an  artificial  radiation  belt  cannot 
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significantly  distort  the  geomagnetic  field.  Moreover,  the  stable -trapping 
limit  exceeds  the  peak  radiation  intensity  for  a one-megaton  belt  by  a 
factor  S 40  at  all  values  of  L . Thus,  one  would  expect  the  Kennel- 
Petschek  limit  to  operate  initially  in  the  case  of  a single  one-megaton 
detonation  if  injection  were  confined  to  a longitude  interval  ^ 10°  in  width. 
The  spectrum  of  fission  electrons,  however,  differs  considerably  from 
the  spectrum  used  by  Kennel  and  Petschek  (1966)  in  calculating  the  limiting 
value  ~7  X 10^^L"^cm~^sec  ^ on  the  integral  omnidirectional  flux  • 

A calculation  of  based  on  (14)  would  thus  be  very  much  in  order, 
before  firm  conclusions  can  be  drawn  concerning  the  relevance  of  the 
Kennel-Petschek  mechanism  to  artificial  radiation  belts. 

The  "Teak"  and  "Orange"  events  of  August  1958  surprised  many 
observers  by  injecting  many  energetic  electrons  onto  quasi-adiabatic 
trajectories.  The  site  conjugate  to  detonation  lay  deep  in  the  dense  atmos- 
phere (in  the  South  Pacific  region),  and  one  might  thus  have  expected 
all  but  a few  of  the  fission  electrons  to  be  in  the  loss  cone.  Davidson 
(1973)  luis  recently  proposed  another  interesting  plasma  instability  to 
account  for  the  discrepancy.  He  notes  that  the  fission-beta  electrons 
from  the  nuclear  debris  formed  a particle  beam  (through  the  ambienv. 
plasma)  directed  toward  the  southern  hemisphere.  The  beam-plasma 
interaction  should  have  generated  electromagnetic  cyclotron  waves, 
according  to  Davidson  (1973),  and  the  resulting  pitch-angle  diffusion 
should  have  enabled  a large  fraction  of  the  fission-beta  electrons  to 
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escape  from  the  loss  cone  (onto  quasi-  adiabatic  trajectories)  before 
reaching  their  southern  mirror  points.  Productive  research  thus  con- 
tinues on  the  topic  of  artificial  radiation  belts,  many  years  after  the  last 
such  belt  ceased  to  exist  (circa  1966). 


6.  EMPIRICAL  STUDIES  OF  PARTICLE  DIFFUSION 


The  purpose  of  analyzing  radiation-belt  observations  is  to  sum- 
marize the  behavior  of  the  particle  distributions  concisely  within  the 
framework  of  known  physical  laws.  This  normally  means  to  specify  the 
numerical  values  of  the  transport  (Fokker -Planck)  coefficients  to  which 
the  particles  have  been  subjected.  One  can,  of  course,  adopt  a "brute- 
force"  approach  and  insert  ail  conceivable  combinations  of  D.  . , D , 
etc.  , into  a numerical  program  that  theoretically  predicts  the  evolution 
of  the  distribution  function  £ (M,  J,9;  t).  From  the  output  one  can  select 
the  prediction  that  best  agrees  with  the  data,  and  thereby  identify  the 
optimal  set  of  transport  coefficients,  i.  e.  , the  set  used  in  generating 
the  best-fitting  prediction. 

Reviewed  here  are  approaches  that  seem  somewhat  more  imagi- 
native, in  that  the  optimum  trial  values  for  D.  . and/or  D are 
extracted  from  the  observational  data  directly.  One  such  approach,  the 
variational  method,  entails  a linear  least- squares  deterpninat,.on  of  that 
D.  . which  makes  D (a  function  of  L)  deviate  minimally  from  a 


-54- 


constant  in  time.  Another,  the  spatial -quadrature  method,  treats  (13) 
as  a first-order  differential  equation  for  In  both  methods,  the 

spatial  and  temporal  derivatives  of  f are  treated  as  given  by  the  obser- 
vational data. 

Having  obtained  trial  values  of  and/or  by  such 

empirical  methods,  one  v/ould  be  prudent  to  test  those  values  in  the  usual 
way,  by  inserting  the  trial  transport  coefficients  into  the  diffusion  equa- 
tion and  seeking  to  predict  the  observed  temporal  evolution  or  stationary 
form  of  f.  However,  the  trial-and-error  aspects  of  the  more  traditional 
approach  are  thus  averted. 

The  data  shown  in  Figure  18  were  extracted  from  satellite  obser- 
vations of  integral  fluxes  in  two  energy  channels  (E>0.  5 MeV  and  E >1.  9 
MeV)  by  assuming  a power-law  energy  spectrum  (Lanzerotti  etal.  , 1970). 
A major  magnetic  storm  began  on  17  December  1962.  This  led  to  a 
redistribution  of  which  continued  for  about  three  weeks,  until 

the  quiet  pre- storm  configuration  was  approximately  restored.  The 
redistribution  was  presumably  accomplished  by  a superposition  of  radial 
diffusion  and  pitch-angle  diffusion,  as  described  by  the  coefficients 
and  respectively.  The  magnitude  of  perhaps  also  of 

D^)  was  presumably  enhanced  during  the  disturbed  period  (e.  g.  , 17-21 
December). 

Lanzerotti  et  al.  (1970)  employed  a variational  method  to  extract 
D,  . and  D from  the  full  set  of  data  described  above.  They  assumed 

Xj  Ij  XX 
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D.  , = D L for  each  value  of  M,  with  D independent  of  L and  time, 
n n 

Defining  F 5 In  f and 


X (L.t)  ■ D 
n ' n 


3F  ^ a^F 

+ y 

3L 


(15) 


they  proceeded  (by  a straightforward  algebraic  operation)  to  minimize 
the  function 

^2  ^2 

s'*-’  j -<V^]  *'*1-  <“> 

Li  tj 

where  g(L)  denotes  an  assigned  weighting  function  and  denotes  the 

temporal  mean  value  cf  X^(L,  t)  at  fixed  L.  The  limits  of  integration  were 
chosen  as  follows:  Lj  = 3.4,  L2  = 4.8,  tj  = 22  December  1962  and  t2  = 

10  January  1963.  The  derivatives  appearing  in  (15)  were  evaluated 
numerically  (using  symmetric  first  differences)  from  the  daily-median 
particle  data. 

By  requiring  dG^/dD^  = 0,  Lanzerotti  et  al.  (1970)  obtained 
D »6  X 10”^®  X day' ^ for  each  value  of  M considered.  The  cor - 

responding  lifetimes  are  shown  in  Figure  19.  The  most 

reasonable  (theoretical)  choice  for  n is  n = 10,  which  would  imply 
D.  T * ^ XIO"^^  day'^.  The  decay  rates  1/t  correspond  to  pitch- 
angle  diffusion.  The  variation  of  t with  L at  fixed  M corresponds  in 
part  to  the  variation  of  energy  with  L at  fixed  M (cf.  Figure  10,  in 
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= 300  MeV/ gauss 


Optimal  electron  lifetimes  obta 
variational  analysis  of  post-sto 
January  1963)  on  outer -zone  eli 


which  T is  an  increasing  function  of  E,  as  well  as  a decreasing  function 
of  L). 

The  functions  "'■(L)  thus  obtained  were  tested  by 

Lanzerotti  et  al.  (1971),  who  integrated  (13)  to  obtain  f as  a function  of 
time  at  each  of  six  L values.  The  initial  configuration  was  established 
by  the  particle  data  at  f (22  December  1962),  and  the  time -dependent 
boundary  conditions  (at  L = 3.  4 and  L = 4.  8)  were  similarly  determined 
from  the  observational  data  at  those  L values.  The  predicted  evolution 
of  f (L,t)  between  L = 3.  4 and  L = 4.  8 was  given  by  a numerical 
integration  of  (13).  The  results,  shown  in  Figure  20  (left  panel),  serve 
to  vindicate  the  values  of  and  t(L)  obtained  for  this  time  interval 

by  means  of  the  variational  method.  The  dashed  curves  represent  the 
prediction;  the  data  points  (connected  by  solid  line  segments)  represent 
the  observation.  A different  set  of  transport  coefficients,  as  proposed 
by  Newkirk  and  Walt  (1968b),  is  similarly  tested  in  the  right  panel  of 
Figure  20.  The  larger  value  of  (=  5 XlO  ^ day”^)  is  plainly 

excessive  for  this  time  interval,  viewed  as  a whole. 

If,  however,  one  considers  the  time  interval  20-31  December 
1962,  the  larger  value  of  » 5 X 10  ^ day  ^)  turns  out  to 

be  preferable.  This  was  the  interval  originally  selected  for  analysis  by 
Newkirk  and  Walt  (1968b).  The  comparisons  between  prediction  and 
observation  are  shown  in  Figure  21,  in  which  20  December  is  used  as 
the  starting  point  (vs.  22  December  in  Figure  20)  for  numerical  integration. 
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numerical  integration  with  time -dependent  boundary  conditions  imposed  by 
the  data  at  L = 3,  4 and  L = 4.  8 (Schulz  and  Lanzerotti,  1974;  based  on 
personal  communication  from  C.  G,  Maclennan,  197Z). 


Figures  20  and  21  thus  suggest  a temporally  decreasing  magnitude 
for  (perhaps  5 XlO”^  day"^  through  25  December  1962,  and 

6 X 10’^*^  day"^  thereafter).  This  should  not  be  surprising,  since 
geomagnetic  conditions  (e.  g. , as  measured  by  the  index  K^)  were  undoubtedly 
more  disturbed  during  the  period  17-25  December  1962  than  during  the 
subsequent  period  between  storms.  It  is  noteworthy  in  this  context  that. 


according  to  Mozer  (1971),  the  spectrum  of  electrostatic  impulses  at 

balloon  altitudes  is  proportional  to  exp  (0.4  ^p)*  electrostatic  impulses 

detected  at  balloon  altitudes  are  indeed  of  magnetospheric  origin,  they 

(C) 

would  thus  lead  to  a radial-diffusion  coefficient  D£^^  that  is  also  propor- 
tional to  exp  (0.  4 Kp).  Lanzerotti  and  Morgan  (1974)  have  made  a similar 
analysis  of  magnetic -disturbance  spectra  observed  on  the  ground,  and 
have  concluded  (within  the  framework  of  the  three-term  Mead  field  model) 
that  magnetic  impulses  produce  a radial-diffusion  coefficient 
exp  (2.  3 Kp)  day"^  for  particles  mirroring  at  the  equator.  Results 
compiled  by  Williams  et  al.  (1968)  suggest  that  D^  varies  much  more 
weakly  than  D.  ^ with  K , and  that  any  such  variation  of  D is  con- 

Xj  P XX 

fined  to  the  "slot"  region  (2.  5 ^ L S 3.  5). 

In  a totally  separate  empirical  analysis,  based  on  the  method  of 
spatial  quadrature,  Farley  (1969)  extracted  a radial-diffusion  coefficient 
D.  , ~ 6 XIO  ^ (1.2/L)®^*^^  day  ^ from  compilations  of  data  (Paulikas 
cl  al.  , 1967;  Imhof  ct  al.  , 1967)  on  inner-zone  electrons.  The  apparent 
decay  rate  (-  9F/8t)  during  the  three  years  following  Starfish  (see  'fable  1, 
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above)  was  considerably  smaller  at  L S 1.  2 than  one  would  have  pre- 
dicted from  the  known  effects  of  Coulomb  collisions.  The  discrepancy  is 
illustrated  in  Figure  22  (right  panel).  Some  process  must  have  partially 
counteracted  the  effects  of  Coulomb  collisions.  Using  the  intensity  profile 
shown  in  Figure  22  (left  panel),  Farley  (1969'  properly  neglected  wave- 
particle  interactions  (i.  e.  , set  = ® ) and  integrated  (13)  for 
Since  (13)  is  a first-order  differential  equation  for  there  is  one 

arbitrary  constant  (taken  by  Farley  to  be  the  value  of  at  L = 1.15). 

The  results,  for  various  choices  of  the  arbitrary  constant,  are  shown  in 
Figure  23  (smooth  curves). 

The  most  acceptable  curve  in  Figure  23  is  the  one  t!.at  corresponds 
-5  -1 

to  = 3 X 10  day  at  L = 1.15;  this  curve  is  rou-jhly  approximated 

by  the  formula  Dj^j^abXlO  (1.2/L>)  day  . The  strong’y  inverse 
Variation  of  with  L had  been  anticipated  by  Newkirk  and  Walt 

(1968a),  whose  earlier  analysis  of  the  same  data  (by  a less  reliable 
method)  had  yielded  the  "staircase"  function  shown  in  Figure  23.  It  is 
difficult  to  imagine  any  physical  process  except  atmospheric  scattering 
that  would  vary  so  strongly  with  L.  Thi^  consideration  led  Roederer 
et  al.  (1973)  to  propose  shell  splitting  by  internal  geomagnetic  multipoles 
(cf.  Figure  8)  as  the  agent  responsible  for  translating  atmospheric  pitch- 
angle  diffusion  into  radial  diffusion  at  very  low  L values.  The  impli- 
cations of  this  suggestion  remain  to  be  worked  out  in  quantitative  detail. 
Moreover,  the  testing  of  this  hypothesis  will  require  the  curves  in 
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Fig,  23, 


Radial  diffusion  coefficients  obtained  from  data  in 
Figure  22,  assuming  fission  spectrum  to  obtain 
f(Mj  J,0)  at  constant  M and  J.  Staircase  function 
(Newkirk  and  Walt,  1968a)  follows  from  self-incon- 
sistent quadrature.  Smooth  curves  are  self-con- 
sistent, but  require  assignment  of  arbitrary  values 
(filled  circles)  to  at  L = 1,15  (Farley,  1969). 
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Figure  23  to  be  replaced  by  solutions  of  the  alternative  equation 


1 _d__ 

l7f  dL 


9F  ^ _1_ 

3t  T 


(17) 


with  the  flux  profile  evaluated  at  fixed  energy,  as  in  Figure  22a.  Newkirk 
and  Walt  (1968a)  had  postulated  a radial-diffusion  process  conserving  M 
and  J,  and  so  had  used  (14)  to  treinsform  the  fltix  profile  in  Figure  22a 
to  a distribution  function  at  constant  M and  J prior  to  numerical  dif- 
ferentiation. One  does  not,  however,  expect  to  find  serious  discrepancies 
(either  in  magnitude  or  in  functional  form)  between  the  solutions  of  (17) 
and  those  of  (13)  for  O.  . . 


7.  INNER- ZONE  PROTONS 

A major  development  in  radiation- belt  theory  (Farley  et  al.  , 1970) 
was  the  inclusion  of  radial-diffusion  effects  in  the  analysis  of  high-energy 
inner-zone  protons.  It  seems  surprising  in  retrospect  that  previous 
authors  (including  Dragt,  1971)  had  chosen  to  neglect  so  pervasive  a 
process  as  radial  diffusion  in  attempting  to  predict  the  spatial  and  spectral 
structure  of  the  inner  proton  belt.  Before  1970,  one  had  vainly  hoped  to 
vmderstand  this  structure  as  a dynamical  balance  beUveen  Coulomb  drag 
(ionization  loss)  and  a source  known  as  GRAND  (the  beta  decay  of  albedo 
neutrons  ejected  from  the  atmosphere  by  incident  cosmic  rays). 
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Farley  et  al.  (1970)  showed  that  the  major  deficiencies  of  such  a 
treatment  are  corrected  by  the  inclusion  of  radial  diffusion  and  the  appro- 
priate boundary  conditions.  The  improvement  is  illustrated  dramatically 
in  Figure  24,  from  a slightly  more  refined  calculation  by  Farley  and  Walt 
(1971).  The  basic  form  of  the  Fokkef -Planck  equation  for  protons, 
neglecting  pitch-angle  diffusion,  is 


where  the  subscript  v designates  a nonstochastic  (frictional)  change  in 
the  invariant  action  integral  Jj^,  the  symbol  S denotes  the  source  term 
derived  from  GRAND,  and  the  symbol  denotes  the  particle  lifetime 
against  charge  exchange  and  nuclear  collisions  in  the  atmosphere.  Charge 
exchange  is  completely  negligible  at  the  high  energies  of  interest  here, 
but  can  be  very  important  at  ring-current  energies  (see  Section  8). 

It  is  convenient  to  transform  (18)  from  the  ' old"  variables  (M,  J,$) 
to  the  "new"  variables  (M,K,L),  where  K^sJ^/bmgM  and  LsZna^  | g|^/® 
The  Jacobian  of  this  transformation  has  the  absolute  vp,.'ue 

|a(M,J,*)/8(M,K,L)|  = (8mgM)^^^  (2ira^/L^)  |g°|,  (19) 


and  the  quantity  K is  conserved  by  the  operative  dynamical  processes. 


Dn-1Cr«L’0day;> 


Planck  equation  with  Dlt  = 0 and  S defined  by  the  Lingenfelter  ( 1963)  spectr 
Solid  curves  (b)  are  steady-state  solutions  of  the  same  Fokker -Planck  equa- 
tion, with  the  same  source  term  S,  but  with  Dll  = 10'°L,*0(jay-l  and  boundary 
conditions  imposed  by  the  data  at  L = 1.1  (where  f = 0)  and  L = 1. 7 (Farley  and 


where  (cf.  Farley  and  Walt,  1971;  White,  1973) 


(dM/dt)^,  = (Y/B^)  (oE/dt)^ 


= (4irqp2  (mQ/2B^)^/2<r 


(21a) 


= <Ng[Y^  - 1 - Y^  In  (Xjjm^v/1i)]> 


+ E <Ni>  - 1 - Y'^  In  [2m^c^  (Y^  - l)/lj|  . 


(21b) 


Here  the  angle  brackets  denote  an  average  over  the  particle  drift  period; 
the  subscript  i denotes  the  atmospheric  constituent  having  density 
(particles  per  unit  volume),  nuclear  charge  number  Z^,  and  mean  excita- 
tion energy  I.  per  atomic  electron;  and  the  symbol  m denotes  the  rest 

X c 

mass  of  an  electron.  The  plasmaspheric  Debye  length  is  equal  to 
2 1/2 

(*fT^/4TrN^q^  ) , where  K is  Boltzmauin's  constant.  The  charge  of  a 

proton  is  denoted  q^. 
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Since  f = J^/Zm^MB,  the  source  term  S must  be  given  (Dragt 
et  al.  , 1966)  by 


S = (1/p^YtJ 


^ j"  • d / 

J HN 


(22) 


where  t (»1013  sec)  is  the  mean  life  of  a neutron  at  rest,  j”  is  the 

n ««» 

differential  neutron  flux  per  unit  solid  angle  in  velocity  space,  and  d / 
is  the  element  of  length  along  a proton's  spiral  path  between  mirror  points. 
Since  inner-zone  drift  shells  are  offset  from  the  geocenter,  the  quantity 
S given  by  (22)  can  vary  with  azimuth.  The  source  term  in  (18)  and  (20) 
should  ideally  be  written  ^SX  so  as  to  indicate  the  drift  average. 

In  obtaining  Figure  24,  Farley  and  Walt  (1971)  set  = m and  solved 
for  the  steady  state  (8f/9t  = 0).  With  ~ equation  (20)  reduces  to 

a first-order  differential  equation  in  M,  having  the  boundary  condition 
that  f = 0 at  M = oo  (actually  imposed  at  M = 4GeV/gauss  for  numerical 
convenience).  The  predicted  distribution  function  (dashed  curves)  is  thus 
deficient  at  low  M (e.  g.  , M = 200  MeV/gauss),  excessive  at  high  M 


(e.  g.  , M = 3 GeV/gauss),  and  generally  in  poor  agreement  with  the  data. 

With  ^ day'^,  one  is  free  to  impose  suitable  boundary  con- 

ditions in  L,  and  the  resulting  predictions  (solid  curves)  for  J^/MB  at 
L < 1.  7 agree  remarkably  well  with  the  observational  data. 

Farley  et  al.  (1972)  added  the  refinement  of  including  the  geomagnetic 
secular  variation  ad  hoc  in  the  basic  equations.  Since  the  secular  variation 
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f 


in  fact  conserves  the  adiabatic  invariants  of  charged-particle  motion,  its 
inclusion  as  a term  in  (20)  can  only  simulate  the  phrnomenon  approxi- 
mately. However,  if  one  views  the  field  configuration  as  fixed  (frozen  in 
time),  then  the  secular  contraction  of  dipolar  drift  shells  can  be  simulated 
by  a radial  convection  of  particles,  an  increase  in  K (to  preserve  the 
equatorial  pitch  angle),  and  a decrease  in  M (to  characterize  the  net 
energy  gain  correctly).  One  must  specify 


(M/M)„^  = (i;/L(„^  = - 2 (k/K)^^^  = gf/g» 


(23) 


for  the  purpose  of  this  simulation.  An  equation  of  the  form 


9f 

— + M 

at 


-l/2_f_ 

8M 


M 


l/2/**^\  1 ..  0 , 0.  , 2 

lirAr  - 


r f 


aL  L L J 


..0/,  0.  3 

- (gi  /2g,  ) — 

* ^ aK 


Kf 


X /•  0/  w-1/2  9 

+ (gi  /gi)  M 

aM 


= L"  — 


1 


df 


3L  LL^  dh 


+ s - — 


(24) 


thus  emerges.  Farley  et  al.  (1972)  obtained  quasi-static  solutions  of 

0 0 

(24)  in  the  limit  K = 0,  taking  the  realistic  case  gj  /gj  = - 1933  yr  and 
(for  comparison)  the  unrealistic  case  g^/g|^  = 0.  The  outcome  is 
shown  in  Figure  25.  The  effect  of  the  geomagnetic  secular  variation  is 
to  enhance  the  predicted  proton  flux  at  each  M,  thus  improving  the 
agreement  with  observation  (notably  at  M = 0.  8 GeV/gauss). 
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Fig,  25.  Inner -zone  proton  distribution  function  (XZnrjQ)  for  J = 0 
and  selected  values  of  M,  based  on  data  points  from 
Figure  24  and  numerical  integration  ( Farley  et  al.  , 
1972).  Solid  curves  are  steady-state  solutions  of  the 
Fokker-Planck  equation  (neglecting  the  geomagnetj^c 
secular  variation)  with - 2,1  y lO'^L^®  X jmQC^/gjM 
and  S defined  by  the  Lingenfelter  ( 1963)  spectrum. 
Dashed  curves  are  quasi-static  solutions  of  the  Fokker- 
Planck  equation  augmented  to  simulate  secular  contrac- 
tion of  the  geomagnetic  dipole  moment  (Farley  et  al,  , 


The  foregoing  results  were  obtained  by  using  an  albedo-neutron 
spectrum  predicted  by  Lingenfelter  (1963)  from  a theoretical  extrapolation 
of  observations  made  at  energies  E<10  MeV.  This  extrapolation  follows 
the  dashed  line  in  Figure  26.  The  subscript  2tt  denotes  angular  integra- 
tion over  the  upward  hemisphere  in  velocity  space.  For  this  presentation 
the  extrapolated  spectrum  of  Lingenfelter  (1963)  has  been  evaluated  at  a 
magnetic  colatitude  6 = 50°.  Until  September  1971,  there  had  been  no 
direct  measurement  of  the  albedo -neutron  flux  in  the  energy  range  E = 
10-100  MeV,  i.  e.  , in  the  range  of  relevance  to  the  majority  of  inner-zone 
proton  observations.  This  deficiency  was  rectified  by  Preszler  et  al. 
(1972),  whose  balloon-borne  neution  detectors  yielded  the  data  points 
shown  in  Figure  26.  It  is  noteworthy  that  the  neutron  fluxes  observed  by 
Preszler  et  al.  (1972)  far  exceeded  the  extrapolation  of  Lingenfelter 's 
prediction.  However,  more  recent  theoretical  (Monte  Carlo)  calculations 
by  Merker  (1972)  and  Armstrong  et  al.  (1973)  have  yielded  results  that 
agree  with  the  measurements  reported  by  Preszler  et  al.  (1972,  1974). 

Fortified  by  these  new  theoretical  and  observational  data,  Claflin 
and  White  (1974)  used  the  CRAND-proton  source  derived  therefrom  in  the 


radial-diffusion  code  that  had  been  developed  by  Farley  et  al.  (1972). 


Clafli.i  and  White  (1974)  included  nuclear  collisions  (the  term  ^/‘q)  in 
their  analysis  and  expressed  the  radial -diffusion  coefficient  as  the  sum 


of  two  terms:  = ^LL  ^ ^LL^  ’ from  electrostatic  impulses 

and  magnetic  impulses,  respectively.  Taking  = 5X10”^^ 


NEUTRON  ENERGY, 


Fig,  26.  Data  points  specify  the  omnidirectional  flux  of  albedo 
neutrons  at  magnetic  latitude  40°N,  altitude  3 6 km, 
as  measured  by  Preszler  et  al.  (1972).  The  dashed 
line  is  an  extrapolation  of  the  spectrum  obtained  by 
Lingenfelter  (1963)  at  energies  E < 10  MeV. 


(M^/M)^  day  ^ and  = B X 10  day  \ where  = 1 GeV/gauss, 

they  obtained  the  theoretical  (dashed  and  solid)  curves  shown  in  Figure  27. 
The  solid  curves  correspond  to  the  atmospheric  model  used  earlier  by 
Farley  and  Walt  (1971).  The  dashed  curves  correspond  to  the  same  neutral- 
atmospheric  model,  but  require  a plasmaspheric  electron  density  five  times 
that  used  by  Farley  and  Walt  (1971).  Agreement  with  observation  is  thereby 
improved,  but  observational  confirmation  of  the  required  plasmaspheric 
electron  density  is  lacking.  One  must  therefore  regard  the  corresponding 
determination  of  and  tentative,  although  the  two-term 

representation  of  in  such  an  analysis  is  definitely  a ' forward. 

A more  detailed  review  of  progress  on  the  problem  of  inner-zone  protons 
is  given  by  White  (1973). 

A word  of  caution  seems  appropriate  at  this  point.  The  use  of 
(24)  to  simulate  the  behavior  cf  inner-zone  protons  is  an  expedient  that 
remains  to  be  justified  in  fundamental  terms.  It  is  true  that  one  can  obtain 
quasi-static  solutions  cf  (24)  by  inc- jrting  transport  coefficients  and  sources 
that  have  been  averaged  over  the  solar  cycle.  One  can  further  obtain 
quasi-periodic  solutions  of  (24)  by  allowing  the  transport  coefficients  and 
sources  to  vary  over  the  solar  cycle  (cf.  Dragt,  1971).  However,  there 
is  no  assurance  that  such  mathematical  solutions  correspond  to  geo- 
physical reality,  even  in  principle,  because  there  is  no  assurance  that 
the  geomagnetic  secular  variation  can  be  treated  properly  in  the  quasi- 
static approximation. 
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Diffusion 

Coefficient 


It  seems  probable,  at  least  for  M >1  GeV/gauss,  that  some  proton 
residence  times  in  the  inner  zone  are  comparable  to  the  8000-year  period 
of  the  geomagnetic  cycle.  If  so,  the  present  configuration  of  the  inner 
proton  belt  would  be  the  product  of  its  past  evolution,  an  evolution  depen-  ? 

dent  upon  the  whole  history  of  the  earth's  field.  (Most  serious  investigators 
reject  the  contrary  hypothesis,  that  the  inner  proton  belt  has  remained 
immutable  since  its  sudden  creation  in  the  distant  past.  ) 

The  proper  evolutionary  treatment  of  inner-zone  protons  follows 
naturally  from  (18).  One  must  attempt  to  model  the  transport  coefficients, 
source  term,  and  boundary  conditions  (in  $)  as  functions  of  time.  All 
consequences  of  the  secular  variation  must  enter  at  this  level,  and  not  as 
diiatiViCt  terms  in  the  Fokker-Planck  equation.  The  distribution  function 
f(M,  J,4;  t)  must  be  stored  as  a function  of  the  adiabatic  invari.jiits  (M,  J,®) 
and  the  time  ( t ).  One  may  presume  that  f(M,  J,$:  ‘)  is  quasi-periodic  on 
the  time  scale  of  8000  years.  Thus,  if  t = 0 deontes  the  present,  one 
may  set  f(M,J,4;  -8000  yr)  = f(M,J,$;  0)  as  an  initial  condition  based  on 
the  observational  data.  One  next  integrates  forward  in  time  (from  t = 

-8000  yr  to  t = 0).  The  chosen  normalization  of  is  thereby  vindi- 

cated  if  the  solution  for  f(M,  J,4i;  0)  matches  the  observational  data  at 
the  present  epoch,  and  the  hist<jrical  description  of  the  inner  proton  belt 
is  thus  completed.  (One  may  instead  integrate  from  the  initial  condition 
that  f = 0 at  the  time  of  the  last  geomagnetic  field  reversal;  see  Cox, 


8.  ALPHA-PROTON  RATIO 


Minor  ionic  constituents  of  the  geomagnetically  trapped  radiation 
are  becoming  recognized  as  a valuable  probe  for  the  dynamics  of  radia- 
tion belts  in  general.  To  the  extent  that  such  minor  ions  (e.  g.  , He^, 

He^^,  O^^)  respond  differently  from  protons  (or  from  each 

other)  to  the  fluctuating  fields  responsible  for  various  transport  processes, 
the  relative  abundances  of  such  ions  can  serve  to  indicate  the  relative 
importance  of  th  i various  transport  mechanisms.  To  the  extent  that  their 
relaiive  abundances  differ  from  the  abundances  inherent  in  various  plau- 
sible radiation  sources,  these  relative  abundances  can  help  one  to  assess 
the  relative  importance  of  the  various  sources. 

It  has  become  conventional,  for  valid  theoretical  reasons,  to 
compare  the  intensities  (fluxes)  of  various  radiation-belt  ions  having  the 
same  velocity,  i.  e.  , the  same  amount  of  energy  per  nucleon  and  (p,*ef- 
erably)  the  same  pitch  angle,  at  a given  location  in  space.  Moreover, 
since  helium  is  the  most  prevalent  minor  ion  in  the  radiation  belts,  the 
alpha/proton  (o/p)  intensity  ratio  is  the  ratio  that  has  been  studied  most 
thoroughly.  (Since  it  would  be  very  difficult  to  distinguish  between  He^ 
tnd  with  instruments  presently  avai  la  bit-  on  spacecraft,  the 

observed  "alpha-p.'irticlc"  intensity  is  understood  to  include  the  contri- 
butions from  both  charge  states.  ) Typical  values  of  the  o/p  ratio 
observed  off-equator  at  L~3  have  been  compiled  by  Krimigis  (1970) 
and  are  shown  in  Figure  28.  The  ratio  is  somewhat  enhanced  (as  shown) 
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during  periods  of  geomagnetic  disturbance  (e.  g.  , November  1968),  but 

_4 

is  limited  even  under  these  conditions  to  values  ~1-30X10  for  E/A  ^ 0.  3 
MeV /nucleon. 

The  magnetospheric  rr/p  ratio  (as  observed  off-equator  at  L~3) 

_2 

is  thus  apprecirbly  smaller  than  the  a/p  ratio  in  the  solar  wind  (~4  X 10~  ), 
which  is  the  most  plausible  source  for  magnetospheric  hydrogen  and  helium 
ions,  Cornwali(1972a)  has  successfully  accounted  for  this  superficial  dis- 
crepancy by  attributing  magnetospheric  radial  diffusion  (off-equator)  to 

(e) 

electrostatic  impulses.  The  resulting  diffusion  coefficient  is  pro- 

2 

portional  to  (Z/A)  at  fixed  energy /nucleon,  among  particles  having  drift 

periods  $ 20  min  (the  postulated  decay  time  of  an  electrostatic  impulse). 

(e)  + 

The  magnitude  of  for  energetic  H (Z  = 1,  A = 1)  thus  exceeds  the 

magnitude  of  for  He^^  (Z  = 2,  A = 4)  and  He^  (Z  = 1,  A = 4)  by 

factors  of  four  and  sixteen,  respectively.  Consequently,  protons  have 
more  efficient  access  than  helium  ions  to  the  inner  magnetosphere  (LS  4) 
icom  an  external  source. 

In  reproducing  the  trend  of  quiet-time  data  shown  in  Figure  28, 
Cornwall  ( 1972a)  neglected  pitch-angle  diffusion,  since  the  observational 
evidence  does  not  yet  require  the  inclusion  of  this  process  for  the  particle 
population  under  study.  He  did,  however,  include  both  ionization  loss 
(Coulomb  drag)  and  charge  exchange,  both  of  which  are  included  in  (20) 
as  far  as  prc  tons  are  concerned.  Ionization  loss  for  helium  ions  is 
described  by  (21)  if  the  correct  ionic  charge  (q^  for  He*^,  Zq^  for  He^^^ 
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is  inserted  in  place  of  the  proton  charge  q^.  Charge  exchange  (’’’q) 
an  especially  important  process  here.  It  leads  not  only  to  the  neutrali- 
zation of  H^,  as  in  (20),  but  also  to  the  alternation  of  helium  ions  between 
the  singly  charged  state  and  the  doubly  charged  state.  To  describe  this 
alternation,  one  requires  a pair  of  coupled  transport  equations  having  the 
form 


■ 1 
7^ 

LL 

dt 

3Lf 

L 
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(25a) 


(25b) 


where  the  total  derivative  (df/dt)  includes  the  Coulomb-drag  term,  as 
on  the  left-hand  side  of  (20).  Distributed  sources  (S)  are  unimportant 
here.  The  charge-exchange  lifetime  x.^  in  (25)  refers  to  the  conversion 
of  He^^  (distribution  function  T)  into  He^'^  (distribution  function  T). 
The  observable  o/p  ratio  is  equal  to  (l6/f  ) (f,  + f^)  at  fixed 

E/A,  where  fsJ^/2mQMB  for  each  species. 

The  foregoing  discussion  refers  to  particles  mirroring  far  off 
the  georr.agnctjc  equator.  The  situation  is  much  different  for  particles 
that  mirror  near  the  equator,  where  Jj^/J^^~10  ^ according  to  the 
initial  observations  by  Fritz  and  Williams  (1973).  Full  equatorial  pitch- 
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angle  distributions  of  inner-zone  alpha  particles  and  protons,  as  obtained 
by  Blake  et  al.  (1973),  are  shown  in  Figure  29.  The  alpha-particle  distri- 


bution is  considerably  narrower  in  pitch  angle  (a)  than  the  proton  distribu- 
tion. This  finding  is  compatible  with  the  previous  observational  results, 
namely  that  oti  the  equator  greatly  exceeds  t'ar  off  the 

equator  along  the  same  field  line. 

Both  pitch-angle  distributions  in  Figure  29  are  sharply  peaked  at 

Q = 9''°.  For  the  lowest  eigenmode  of  pitch-angle  diffusion,  one  might 

2 

have  ei  pected  a distribution  like  sin  a or  sin  a.  It  is  natural,  therefore, 
to  neglect  pitch-angle  diffusion  in  seeking  to  explain  the  observed  distri- 
butions. Nakada  et  al.  (1965)  confronted  the  same  problem  years  ago,  and 
concluded  that  radial  diffusion  from  an  external  source  is  sufficient  to 
produce  a sharply  peaked  pitch-angle  distribution  at  L$3.  The  point  is 
that  particles  from  a common  source  gain  relatively  more  energy  per  unit 
change  in  L at  small  values  of  K than  at  large  values  of  K.  The  dashed 
curves  in  Figure  30  (which  converge  to  a common  value  of  E «9  keV  at 
Lal2)  illustrate  this  effect  for  nonrelativistic  particles.  The  specified 
angles  (90°,  30°,  20°,  10°)  represent  the  equatorial  values  of  o at  L = 7. 
Thus,  particles  contributing  to  a given  energy  channel  at  L~2  (Figure  2S’) 
at  small  a (large  K)  must  have  ari.  en  from  a higher  energy  on  the  source 
spectrum  than  particles  contributing  to  the  same  channel  at  a = 90°  (K  = 0). 

Since  the  source  (plasma-sheet)  spectrum  has  an  £-folding  energy 
of  only  <6  keV  for  protons,  it  is  a steeply  falling  function  of  energy  for 
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Fig,  29.  Equatorial  pitch -angle  distributions  of  alpha 
particles  (left  panel)  and  protons  (right  panel) 
in  the  same  four  energy/ nucleon  passbands, 
as  observed  on  spacecraft  OVl-19  (Blake  et  al 
1973).  Statistical  error  bars  would  be  about 
as  large  as  the  data -point  symbols  themselves 


.6  2 3 4 5 6 7 

L 

Fig,  30.  Empirical  e -folding  energies  (solid  curves)  of 
observed  proton  spectra  (Davis  and  Williamson, 
1963)  at  equatorial  pitch  angles  consistent  with 
conservation  of  M and  J;  expected  variation 
(dashed  curves)  of  energy  with  L for  individual 
protons  having  constant  M and  J,  for  selected 
values  of  energy  and  equatorial  pitch  angle  at 
L = 7 (Nakada  etal. , 1965). 
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E ^ 20  keV.  This  feature  is  translated  by  radial  diffusion  into  a sharply 
peaked  pitch-angle  distribution  at  sufficiently  high  energies  and  low  L 
values.  Nakada  et  al.  (1965)  expressed  the  same  result  somewhat  dif- 
ferently, in  terms  of  the  e^-folding  energies  Eq  of  the  exponential  spectra 
observed  at  various  values  of  K and  L.  Using  the  proton  data  of  Davis 
and  Williamson  (1963),  they  made  the  empirical  observation  that  the 
spectral  index  Eq  varies  with  L in  roughly  the  same  manner  as  the 
kinetic  energy  of  an  individual  particle  would  vary.  This  property  is 
illustrated  by  the  solid  "curves"  in  Figure  30,  where  the  value  of  K is 
uniquely  determined  by  specifying  the  equatorial  value  of  a (90°,  30°, 

20°,  10°)  at  L = 7. 

Under  very  restrictive  conditions,  one  can  prove  a theorem  which 
holds  that  the  £- folding  energy  of  an  exponential  spectrum  of  nonrelativistic 
particles  should  indeed  behave  as  the  energy  of  an  individual  particle.  The 
proof  requires  neglect  of  Coulomb  loss  (dM/dt)^  and  distributed  sources 
(S),  the  assumption  of  a steady  state  (9f/9t  = 0),  and  the  existence  of  an 
M-independent  inner  boundary  (Lj)  at  which  f vanishes.  Then,  if 
T D.  , is  independent  of  M,  the  solution  of  (20)  for  f at  fixed  K factors 
such  that  f(L,  M)  = g(L)  h(M).  Sho  ild  h(M)  be  of  the  form 


t'i 


h(M)=M  exp(-M/MQ), 


(26a) 


the  spectrum  of  nonrelativistic  energies  will  vary  with  L as 
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(26b) 


I 


I 


h(M)  = (B^/E)'  exp(-E/EQ). 

where  E^  = M^B^.  Thus,  the  ^-folding  energy  Eq  varies  with  L as 

would  the  energy  of  an  individual  particle  conserving  M (=Mq)  and  K, 

thereby  tracing  an  L-dependent  mirror-point  field  This  proof 

includes  the  limiting  cases  t = oo  (negligible  charge  exchange),  = ® 

^ u 

(power-law  energy  spectrum  for  J^^),  and  / = 1 (exponential  energy 
spectrum  for  J^).  The  most  doubtful  asumption  (if  < co)  is  the  M- 
independence  of 

It  remains  to  understand  why  the  alpha-particle  distribution  is 
even  narrower  in  pitch  angle  than  the  proton  distribution  (Figure  29). 

One  possibility,  suggested  by  Cornwall (1972a)  in  a different  context,  is 
that  the  "temperature"  of  the  alpha-particle  source  is  less  than  four  times 
the  "temperature"  of  the  proton  source.  Although  E/A  is  the  same  for 
alpha  particles  and  protons  in  the  solar  wind  up  to  the  bow  shock,  it  is 
at  least  conceivable  that  the  two  species  would  exchange  energy  (partially 
thermalize)  in  the  turbulent  magnetosheath.  Further  thermalization  in 
the  pl.isma  sheet  hardly  seems  avoidable,  in  view  of  the  turbulence 
required  there  to  dissipate  the  neutral-sheet  current  and  thereby  main- 
tain the  dawn-dusk  magnetospheric  electric  field. 

Another  possibility  is  that  radial  diffusion  by  magnetic  impulses 
can  affect  the  pitch-angle  distribution  of  alpha  particles  (more  so  than 

for  protons).  The  rationale  for  this  thought  is  illustrated  schematically 

(m) 

in  Figure  31.  The  known  functional  form  of  , as  given  by  Falthammar 
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Fig,  31,  Schematic  representation  of  pitch-angle  and  charge - 
state  dependence  of  radial-diffusion  coefficient  re- 
sulting from  step-like  electrostatic  impulses  (solid 
curves)  or  magnetic  impulses  (dashed  curve).  The 
relative  displacement  of  right  and  left  ordinates 
remains  to  be  determined  empirically  (Fennell 
et  al,  , 1974)  for  given  E/A  and  L, 
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(1968),  strongly  favors  the  access  of  particles  having  equatorial  pitch 

angles  o-q  near  90°.  The  magnituJ.'  of  is  independent  of  particle 

species,  charge  state,  and  energy.  Conversely,  the  known  functional 

(c) 

form  of  D£^j^  slightly  favors  the  access  of  particles  having  small  values 

(e) 

of  sin  Oq.  At  fixed  energy/nucleon,  the  magnitude  of  varies  as 

2 Iq\ 

(Z/A)  for  drift  periods  £ 20  mir..  Since  also  var’es  as 

2 2 

(Y/p  ) for  fixed  Z and  A,  there  should  be  a range  of  E/A  values  for 

which  the  functions  overlap  as  illustrated  (Figure  31). 

This  condition  would  tend  to  make  the  alpha-particle  distribution  narrower 

in  pitch  angle  than  the  proton  distribution. 

Lest  there  be  any  doubt  that  inner-zone  alpha  particles  in  fact 

have  an  external  source,  a near-equatorial  profile  of  2mQf  at  fixed  K 

2 

and  M is  shown  in  Figure  32.  (A  convenient  unit  for  measuring  K and 

$ is  the  giga-weber,  GWb.  For  example,  a particle  having 
2 

at  L~2  yields  K ~3  GWb  and  $ 4 GWb.)  The  diffusion  "current", 

whose  "divergence"  appears  in  (20),  is  equal  to  -Dj^j^(9f /9L)j^  It  is 
directed  inward.  It  is  found  (Blake  et  al.  , 1973)  that  a diffusion  coefficient 
?s  1.  2 X 10'®  day  ^ would  be  required  to  maintain  the  profile 
shown  in  Figure  32  against  Coulomb  loss,  so  as  to  make  9f/9t  vanish. 

Although  it  is  tacitly  assumed  in  several  of  the  paragraphs  above 
that  radiation-belt  particles  are  ultimately  derived  from  the  solar  wind, 
this  assumption  has  not  yet  been  verified  conclusively.  Thus,  Axford 
(1970)  has  suggested  that  perhaps  the  earth's  plasma  sheet  is  populated 
also  by  ionospheric  plasma,  which  is  believed  to  escape  along  open  field 
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lines  in  the  form  of  a "polar  wind"  (e.  g.  , Banks  and  Holzer,  1968). 

Various  tests  have  been  proposed  to  distinguish  between  the  solar-wind 
and  polar-wind  source  hypotheses.  Axford  (1969i  1970)  has  noted  that 
the  two  sources  differ  significantly  with  respect  to  isotopic  and  charge- 
state  abundances  of  ionic  helium.  Mogro-Campero  and  Simpson  (1970) 
have  noted  that  they  differ  with  respect  to  the  ionic  abundances  of  carbon, 
nitrogen,  and  oxygen  (CNO).  Blake  (1973)  has  determined  that  a measure- 
ment of  the  C/O  intensity  ratio  in  the  radiation  belts  would  provide  the 
most  decisive  test,  since  the  two  proposed  sources  differ  by  at  least  five 
orders  of  nvagnitude  in  their  C/O  abundance  ratios.  Observations  reported 
by  Mogro-Campero  (1972)  indicate  that  C/O~l-10  at  L = 4,  for  E/A~  15-30 
MeV/nucleon.  This  result  would  exclude  the  polar  wind  (in  which  C/O  <10 
as  a source  of  such  radiation-belt  ions,  and  favor  the  solar  wind,  in  which 
C/0~1  (Blake,  1973).  However,  the  C/O  intensity  ratio  remains  to  be 
measured  at  the  lower  values  of  E/A  that  would  be  compatible  with  Figure 
28,  i.  e.  , at  energies  which  trace  back  (at  constant  M and  K)  to  the 
nearly  thermal  component  of  the  plasma  sheet. 

Finally  it  is  of  interest  to  comment  on  the  radiation-belt  ions  that 
are  suddenly  lost  through  charge  exchange  at  L ^ 3,  in  particular  on  their 
reappearance  as  a partial  radiation  belt  at  L values  between  1.00  and 
1.15.  Hovestadt  et  al.  (1972b)  have  observed  such  a partial  belt  of  protons 
at  E > 500  keV,  while  Mizera  and  Blake  (1973)  have  extended  the  obser- 
vations down  to  12.  4 keV  (i.  e.  , through  energies  typical  of  the  ring 
current).  The  formation  of  such  partial  radiation  belts  has  been  explained 
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by  Moritz  (1972).  Charge  exchange  converts  an  energetic  ion  (H^,  for 
example)  into  a fast-moving  neutral  atom,  whose  trajectory  is  tangential 
to  the  adiabatic  spiral  path  of  the  ion  at  the  instant  of  neutralization.  Most 
such  rectilinear  trajectories  allow  the  fast  neutral  atom  to  escape  from 
the  magnetosphere,  but  a certain  fraction  are  directed  toward  the  dense 
atmosphere.  In  this  latter  case,  there  arises  a significant  probability 
for  the  dense  atmosphere  (e.  g.  , at  altitudes  ^ 120  km)  to  strip  away  the 
fast  atom's  electron,  thus  recreating  an  energetic  ion.  Ions  thus  formed 
at  low  altitudes  drift  westward  to  form  a partial  radiation  belt,  but  are 
ultimately  degraded  by  Coulomb  loss.  Successive  charge  -exchange  and 
stripping  reactions  within  the  dense  atmosphere  can  lead  to  a novel  form 
of  radial  diffusion.  However,  ions  trapped  at  very  low  L values  (L  ^ 1.  1) 
cannot  complete  a drift  period  around  the  earth.  They  are  doomed  to  lose 
their  energy  at  the  longitude  of  the  South  Atlantic  "anomaly",  if  not  before. 
The  equatorial  pitch-angle  distributions  of  such  recreated  ions  are  sharply 
peaked  about  90°,  because  the  source  protons  at  L 2 3 are  similarly 
anisotropic  (Moritz,  1972;  Mizera  and  Blake,  1973). 

9.  COLD-PLASMA  INJECTION 

Recent  theoretical  and  observational  studies  suggest  rather  strongly 
that  the  plasmapause  (see  Figure  2)  is  a boundary  of  great  importance  in 
the  context  of  magnetospheric  wave-particle  interactions.  This  importance 
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arises  mainly  from  the  fact  ^hat  momentum- space  instabilities  of 
radiation-belt  particles  typically  have  growth  rates  that  are  very  sensitive 
to  the  local  density  of  cold  plasma.  Thus,  the  electromagnetic  ion- 
cyclotron  waves  responsible  for  precipitation  of  ring-current  protons 
(Cocke  and  Cornwall,  1967;  Cornwall  et  al.  , 1970)  and  for  the  electron 
heating  prerequisite  to  the  formation  of  stable  auroral  red  (SAR)  arcs 
(Cole,  1965;  Cornwall  et  al.  , 1971)  in  the  atmosphere  should  be  unstable 
(spontaneously  generated)  only  within  the  plasmasphere.  Ring-current 
observations  analyzed  by  Russell  and  Thorne  (1970)  and  SAR-arc  obser- 
vations analyzed  by  Chappell  et  al.  (1971)  seem  to  confirm  this  prediction. 

The  knowledge  that  cold-plasma  density  should  similarly  affect 
electron-cyclotron  instabilities  (e.  g.  , Kennel  and  Petschek,  1966)  led 
Brice  (1970,  1971)  to  propose  th^^experimental  injection  of  artificial 
plasma  clouds  into  the  magnetosphere  as  a means  of  simulating  sucn 
wave -amplifying  properties  of  the  plasmasphere.  The  effect  of  increasing 
the  cold-plasma  density  (N^)  is  to  decrease  the  phase  velocity  (w/k||)  of 
an  electromagnetic  cyclotron  wave  propagating  in  the  whistler  mode. 

This  decreases  the  electron  energy  required  for  resonance,  in  the  sense 
of  (8),  with  a wave  of  frequency  w/Zir.  Since  the  lower-energy  electrons 
are  typically  the  more  numerous,  the  magnitude  of  the  growth  rate  is 
thereby  enhanced.  Cuperman  and  Landau  (1974)  have  investigated  this 
effect  quite  thoroughly  for  a bi -Maxwellian  plasma  (of  density  Nj^)  in 
which  T^^  = (A  + 1)  T||^  . The  parameter  A ( > 0)  is  known  as  the  anisotropy. 
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If  (=  8TrNj^#iT^®/B^)  is  much  less  than  unity, 

the  growth  rate  Im  w is  well  approximated  by  the  formula 

Im  o)  = -VgNj^  (q/nw)^[Q  + ( w -n)  (A+1]J  (2ir^ 

X exp  [.  (ni^c^/Zn^cu^AfT®)  (.T-w)^],  (27) 

where  n is  the  refractive  index  and  v is  the  group  velocity.  The 

s 

former  is  given  by 

n^  = 1 + (4TrN^q^/m^u)  (fl-u)"^  » 1,  (21) 

-3 

where  f2/2ir  is  the  electron  gyrofrequency  and  w ^ 10  fl . For  a given 

wave  frequency  w /2tt  , the  temporal  growth  rate  Jm  w is  maximized  in 

absolute  value  (Cuperman  and  Landau,  1974)  by  setting  the  argument  of 

the  exponential  equal  to  1/2,  The  more  relevant  spatial  growth  rate 

(Im  cj/v  ) is  maximized  in  absolute  value  by  setting  the  arg-.’ment  of  the 
6 

exponential  equal  to  unity.  This  requires  the  condition 

Nj^/N^  = (fi-w)'^.  (29) 

Since  lmw>0  requires  w /SI  < A/(A+1),  the  optimum  cold-plasma  density 

2 G • 

for  wave  growth  has  the  property  that  N^/N^<  Ap^'  (A+1)  . For  P||  « i, 
the  frequency  corresponding  to  the  maximum  value  of  Im  w/v  lies  just 


slightly  below  w/2tt  = [A/(A+l)j(fl/2ir),  and  so  the  optimum  ratio  Nj^/N^ 


e 2 

is  in  fact  a number  of  order  Ap„  (A+i)  . 


Since  realistic  conditiono  often  invalidate  the  foregoing  assumptions 

(e.  g. , that  « 1 and  that  Nj^«kN^),  Cuperman  and  Landau  (1974)  have 

extended  their  calculation  analytically  to  conditions  far  more  generalthan 

those  admitted  by  (27)-(29).  They  find,  for  example,  that  the  value  of 

Im  w (not  Im  w/v  ) is  maximized  (with  respect  to  both  N and  w)  by  setting 
8 c 

[l  + (N^/Nj^)]P®  A(A  + 1)^  = 2.  Thus,  the  presence  cf  • old  plasma  (N^) 

e 2 

enhances  the  maximum  temporal  growth  rate  (Im  w)  if  p^  A(A+1)  < 2. 

Cuperman  and  Salu  (1974)  have  checked  this  and  other  analytical  results 
nume  rically. 

The  calculations  described  above  invoke  an  energy-independent 
anisotropy  A.  Lucas  and  Brice  (1973)  have  called  attention  to  the  effects 
of  an  energy-dependent  anisotropy.  They  have  pointed  out  that  on  the 
distant  night  side  of  the  magnetosphere  (Figure  7,  left  panel),  the  pitch- 
angle  anisotropy  at  high  particle  energies  is  strongly  unfavorable  for 
wave  growth.  By  reducing  the  particle  energy  required  for  cyclotron 
resonance  at  a given  frequency  w/2ir,  the  addition  of  cold  plasma  in  this 
region  of  space  would  enable  the  wave  to  interact  with  electrons  having 
a more  favorable  anisotropy  (A>0)  for  inducing  wave  growth.  Detailed 
calculations  (Lucas  and  Brice,  1973)  have  supported  this  heuristic  argu- 
ment. Moreover,  the  natural  plasma  in  this  region  of  space  is  so  tenuous 

_3 

(Nc  « N^~l  cm  ) that  even  a modest  experimental  injection  of  cold 

_3 

plasma  (AN^~3  cm  , for  example)  would  constitute  a large  relative 
perturbation  of  the  medium. 
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The  growth-rate  enhancement  of  an  electron-cyclotron  wave  does 
not  depend  on  the  ionic  species  of  the  injected  plasma.  Thus,  barium 
and  cesium  plasmas  are  just  as  effective  (per  unit  density  N^)  as  lithium 
and  hydrogen  plasmas  in  promoting  electron- cyclotron  instabilities.  How- 
ever, the  electron  population  contains  only  a minor  fraction  of  the  total 
magnetospheric  particle  energy.  The  great  majority  of  magnetospheric 
particle  energy  is  concentrated  in  the  population  of  ring -current  protons 
(E  — 10-100  keV).  Phase  velocities  of  electromagnetic  ion-cyclotron 
waves  are  reduced  (and  growth  rates  thereby  enhanced)  by  cold-plasma 
injection  only  at  frequencies  below  the  gyrofrequency  of  the  injected  ion. 
Since  the  natural  value  of  Im  u/v^  attains  its  maximum  at  an  appreciable 
fraction  of  the  critical  frequency  u /2ir  = Anp/2ir(A+l),  efficient  enhance- 
ment of  Max  Im  w/v  by  cold-plasma  injection  requires  the  use  of  a 

S 

fairly  light  ion  (M&rk,  1974).  Moreover,  the  dimensions  of  the  plasma 
cloud  must  exceed  the  wavelength  of  the  alleged  instability.  Technological 
considerations  related  to  this  last  point  led  Cornwall  (1972b)  to  dismiss 
barium-  and  cesium-plasma  injections  in  favor  of  lithium-plasma  injec- 
tion as  a means  of  tapping  the  free  energy  inherent  in  the  proton  ring 
current.  (Still  lighter  ions  such  as  helium  and  hydrogen  cannot  be  pro- 
duced in  sufficient  numbers  with  conventional  energy  sources  on  space- 
craft. Substances  that  are  easily  photo-ionized  by  the  sun,  e.  g.  , 
alkali  metals,  circumvent  this  technological  difficulty.  ) 

The  calculation  of  growth  rates  Y ( s Im  w)  in  the  earth's 
natural  ring  current  (located  beyond  the  plasmasphere)  is  complicated 
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by  the  virhial  absence  of  "old  plasma  and  by  the  fact  that  >«<  1 there. 
Thus,  a treatment  modeled  after  (27)-(29)  is  in- 
appropriate, and  a treatment  based  on  the  full  plasma  dispersion 

equation  is  ultimately  required.  Indeed,  if  the  cold-plasma  dispersion 
relation 


n^  = 1 + (c/v^)^  [l  - (w/fi^)]"^  »1 


(30) 


is  used  in  order  to  compute  the  phase  velocity  w/k  for  k parallel  to 
the  result  is  that  shown  by  the  dashed  curve  in  Figure  33,  where  v^  is 
the  Alfven  speed.  The  temporal  growth  rate  that  would  follow,  assuming 
a sharp  resonance  at  the  value  of  v given  by  (8),  is 


Y = 2N  V (irq/ncw)^  (2Trm  [AO  - (A+l)w] 

p g P * 


X exp  [-  («p-w)^/p„w^np] 


(31) 


if  the  plasma  electrons  are  cold.  This  result  is  shown  as  the  dashed 
curve  in  Figure  34,  Howeve",  the  correct  phase  velocity  Re  w/l^n  and 
growth  rate  Im  w for  a bi-Maxwellian  plasma  are  given  by  solutions 
of  the  dispersion  equation 


c k, 


2 2 
w + (4ttN  q w/n  m ) 
e'l  p p' 


+ (4ttN  q^/m  )|A  - [AO  - (A  + l)w] 
P P P 


x(mp/2ACTPk2)l/2 


(32) 
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Fig.  34.  Normalized  growth  rates  of  electromagnetic  ion-cyclotron 
waves  from  exact  numerical  calculation  (solid  curve)  and 
from  lowest-order  (low-beta)  determination  of  phase  veloc 
ity  (dashed  curve).  The  exact  treatment  increases  the 
phase  velocity  (cf.  Figure  33),  and  thereby  reduces  the 
growth  rate,  at  frequencies  w/2it  < Anp/2ir(A+l),  This 
figure  originally  appeared  with  incorrect  normalization 
in  Cornwall  and  Schulz  (1971).  Normalization  was  cor- 
rected by  Cornwall  and  Schulz  ( 1973). 


where  Reu<n^  and  where  Z(^)  is  the  plasma  dispersion  function  of 
Fried  and  Conte  (1961).  The  (complex)  argument  ^ is  equal  to 
{vn^/ZkT^k^)  (w-Op).  The  resulting  values  of  Re  w/k||V^  and 
Im  w/Op  (where  positive)  are  plotted  as  solid  curves  in  Figures  33  and 
34  (Cornwall  and  Schulz,  1971;  1973;  based  on  personal  communication 
from  T.  Samec,  1971).  The  exact  and  (lowest- order)  approximate  phase 
velocities  (Figure  33)  agree  at  marginal  stability,  i.  e.  , at  w = AJ2p/(A+l), 
because  the  square -bracketed  coefficient  of  Z(^)  in  (32)  vanishes  there. 

The  quantitative  consequences  of  injecting  various  species  of  cold 
plasma  into  such  a proton  ring  current  have  been  investigated  by  Mark 
(1974).  His  results  are  shown  in  Figure  35.  The  phase  velocities  in 
this  case  were  calculated  by  retaining  the  first  three  termc  (n  = 0,  1,  2) 
of  the  asymptotic  expansion 

CO 

CZ(U  - - r(n  + l),  (33) 

n=C  2 

which  holds  for  Im  ^ > 0 (Fried  and  Conte,  1961).  Retention  of  only 
the  first  term  (n  = 0)  leads  to  the  "low-beta”  approximation  shown  in 
Figures  33  and  34.  The  validity  of  Mark's  "second-order”  approxima- 
tion (expansion  of  ^Z  through  n = 2 rather  than  n = w ) can  be  checked 
by  comparing  the  peak  values  of  V/flp  in  Figures  34  and  35  for  the  case 
of  absent  cold  plasma.  One  thus  obtains  Max  Y/Jlp  = 0.  100  by  retaining 
only  one  term  in  (33),  Max  Y/flp  = 0.  035  by  retaining  three  terms 
(Fig..re  35),  and  Max  Y/Op  = 0.044  by  retaining  all  terms  (Figure  34). 
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Fig.  35,  Normalized  growth  rates  of  electromagnetic  ion- 

cyclotron  waves  in  plasmas  consisting  of  hot  protons, 
cold  electrons,  and  (in  two  cases)  additional  cold 
plasma  at  10%  of  the  proton  number  density  (personal 
communication  based  on  Mark,  1974).  Asymptotic 
expansion  of  the  plasma  dispersion  function  (Fried 
and  Conte,  1961)  was  truncated  after  three  terms  in 
this  calculation  (c£.  Figure  34  for  only). 


The  addition  of  cold  plasma  (in  Figure  35,  a 10%  increase  in  N^) 

enhances  the  growth  rate  at  w < as  expected,  and  reduces  the  growth 

rate  at  u > O , where  fl  /2ir  is  the  gyrofrequency  of  the  ionic  additive. 

^ c 

However,  since  Max  V/ft  occurs  at  w ~ 0.  4 i:»  the  absence  of  cold 

P P 

plasma  for  this  situation  ( = A = 1),  there  is  no  technologically  feasible 

ionic  additive  that  can  act  to  enhance  the  existing  peak  in  Y/f2p.  One  must 

hope  instead  to  create  a new  maximum  in  Y/f2  at  u>  < n (by  means  of 

p c ' 

a cold-plasma  admixture  » 10%).  The  natural  maximum  in  Y/fl  would 

P 

occur  in  the  vicinity  of  for  Li^)  only  for  rather  small 

values  of  A and/or  rather  large  values  of  p , i.  e.  , under  conditions 
that  are  already  rather  disturbed  even  in  the  absence  of  cold-plasma 
additives.  Thus,  the  conditions  most  favorable  to  efficient  experimental 
enhancement  of  Y/Jl  are  the  conditions  most  likely  to  preclude  unequivocal 
detection  of  such  enhancement. 

For  this  reason  and  others,  the  idea  of  injecting  cold  plasma  to 
stimulate  electromagnetic  ion-cyclotron  instabilities  has  (in  recent  years) 
fallen  into  some  disfavor  as  an  experimental  project.  Besides  implying 
the  above  objection,  Matk  (1974)  has  argued  that  the  natural  value  of 
Max  Y/v  typically  exceeds  the  reciprocal  plasma  dimension  (~l/La) 

O 

even  beyond  the  plasmasphere.  Thus,  local  wave  growth  would  exceed 
convective  wave  loss,  and  there  would  be  instability  even  in  the  absence 
of  a cold-plasma  additive  (a  point  conceded  by  Cornwall  and  Schulz,  1973). 
Moreover,  Coroniti  et  al.  (1972)  have  shown  that  ring -current  protons 
having  Pu  ^ A ^ 1 beyond  the  plasmapause  should  destabilize  a 
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quasi-electrostatic  wave  mode  (e.  g. , Post  and  Rosenbluth,  1966).  Indeed, 
Burch  (1973)  and  Mizera  (1974)  have  observed  precipitating  ring-current 
protons  in  strong  pitch-angle  diffusion  beyond  the  plasmasp>  ere.  Such 
protons  are  probably  responsible  for  hydrogen  arcs.  The  electrostatic 
loss-cone  mode  is  stabilized  by  an  increase  in  the  cold-plasma  density 
(Coroniti  et  al.  , 1972).  Thus,  the  pl-'smapause  should  serve  as  a boundary 
between  regions  in  which  electromagnetic  and  electrostatic  ion-cyclotron 
instabilities  are  respectively  predominant.  However,  the  experimental 
consequences  of  cold-plasma  injection  are  much  more  subtle  for  ring- 
current  protons  than  one  had  originally  expected  (e.  g.  , Cornwall,  1972b). 

A concise  review  of  present  uncertainties  in  this  evrea  of  investigation,  as 
applied  to  the  natural  environment,  is  given  by  Coroniti  (1973). 

Theoretical  problems  posed  by  the  idea  of  cold-plasma  injection 
include  the  identification  of  a limit  on  stably  trapped  particle  flux,  the 
concept  described  by  Kennel  and  Petschek  (1966).  This  limit  is  the  sup- 
posed maximum  particle  intensity  that  can  remain  trapped  in  the  geo- 
magnetic field  without  provoking  the  spontaneous  excitation  of  large - 
amplitude  electromagnetic  cyclotron  waves.  Its  evaluation  is  based  on 
the  requirement  that 

R exp(2YLa/v  ) < 1,  (34) 

S 

where  La  is  the  effective  path  length  over  which  wave  amplification 
occurs,  and  R is  the  effective  reflection  coefficient  to  which  the  wave 


f 


[ 

\ 

I 


\ 

r 

I 


-102- 


energy  is  subjected  at  the  "end"  of  the  amplification  path.  Stability 


requires  that  the  upper  bound  on  Y/v  , namely  (l/2La)|ln  r|,  be  i.nposed 

S 


at  all  relevant  wave  frequencies  w/2ir. 


If,  following  Kennel  and  Petschek  (1966),  one  imposes  a pitch- 


angle  distribution  (sin  a,  with  A*^l  and  energy  spectrum  (E~”,  with 


n*«4)  at  the  outset,  then  one  derives  a limiting  magnitude 


(E  ) ^7X10  L (Bq/B)  cm  sec 


for  the  integral  omnidirectional  flux  above  the  critical  energy  E . The 


critical  energy  for  evaluating  is  given  by 


E*  = Bq /87tN^(A+1)^A 


(36a) 


for  electrons  and 


E*  = Bq /8TrNp(A+l)A^ 


(36b) 


for  protons  (e.  g.  , Cornwall,  1972).  Thus,  the  addition  of  cold  plasma 
reduces  E . The  value  of  B/Bq  in  (35)  determines  the  location  along 
a field  line  (L),  with  B/Bq  = 1 on  the  equator  and  B/B^  > 1 off  equator. 
If,  following  Brice  and  Lucas  (1971),  one  imposes  only  the  pitch- 


angle  distribution  as  above,  then  a limiting  (nonrelativistic)  spectrum 


of  the  form 


E J4*(E)  - 10^®  (Bq/B)^  cm’^  sec"^ 


w 
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emerges  for  E »E  . Here  the  symbol  J^^(E)  denotes  the  differenti.'.l 

){c  )Jc 

omnidirectional  particle  flux.  The  integral  flux  (E  ) derived  from 
(37}  would  diverge  logarithmically.  There  is  no  possibility  of  reconciling 
(37)  with  (35),  since  different  assumptions  about  the  particle  spectrum 
have  been  invoked  in  deriving  the  respective  equations. 

Both  of  the  foregoing  procedures  have  been  criticized  by  Etcheto 
et  al.  (1973),  who  argue  that  neither  the  energy  spectrum  nor  the  pitch- 
angle  distribution  should  be  specified  a priori  in  a theoretical  calculation 
of  the  other.  Instead,  they  advocate  the  steady- state  solution  of  coupled 
equations  for  the  wave  spectrum  *^^(fa>/2ir)  and  particle  distribution  func- 
tion f (v  , V ) in  the  presence  of  a particle  source  S.  Similarly,  Haerendel 
(1970)  views  the  observable  pitch-angle  distribution  and  energy  spectrum 
as  a consequence  of  the  balance  between  radial  diffusion  and  pitch-angle 
diffusion.  The  main  difficulty  in  implementing  these  approaches  has  been 
that  the  waves  were  assumed  to  propagate  only  parallel  to  ^ Thus,  for 
reasons  outlined  in  Section  3,  the  resulting  pitch-angle  distributions  have 
tended  to  form  a cusp  (A  = + ® ?)  at  Oq  = 90°.  This  feature  makes  it 
impossible  to  define  E , and  (of  course)  disagrees  with  the  observational 
data. 

It  is  evident,  however,  that  the  advocates  of  sslf-consistently 
determining  f(v  , v ) and '6  (w/2tr)  are  on  the  right  track.  A similar 

n 1 1 

procedure  had  been  proposed  by  Cornwall  (1966).  What  remeiins  is  to 
model  the  multitude  of  interactions  more  realistically,  as  Lyons  et  al. 

(1972)  have  done  for  the  case  of  parasitic  pitch-angle  diffusion.  Thus, 
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one  must  include  off- equatorial  resonances  and  oblique  wave  propagation 
(as  well  as  radial  diffusion  of  particles)  in  order  to  obtain  a realistic 
distribution  function.  It  will  be  inevitable  in  a realistic  treatment,  there- 
fore, that  waves  generated  on  a given  L shell  will  interact  later  with 
particles  on  another.  The  various  drift  shells  are  not  dynamically  isolated, 
even  with  respect  to  pitch-angle  diffusion.  Rather,  the  entire  magneto- 
sphere (including  the  ionosphere,  according  to  the  most  thoughtful  investi- 
gators) constitutes  a single  dynamical  entity. 

Pending  a theoretical  solution  of  the  whole  problem,  however, 
there  can  be  no  serious  objection  to  the  use  of  observed  pitch-angle  distri- 
butions and  particle-energy  spectra  in  estimating  the  consequences  of 
cold-plasma  injection.  This  was  the  procedure  adopted  by  Lucas  and 
Brice  (1973),  for  example.  It  is  always  hazardous,  of  course,  to  predict 
the  consequencf.4^  of  a major  environmental  perturbation  before  a full 
understanding  of  the  natural  equilibrium  ha&  evolved.  It  is  therefore  to 
be  hoped  that  future  efforts  in  space  will  be  directed  as  much  toward  an 
understanding  of  the  natural  magnetosphere  as  toward  methods  of  its 
optimal  contamination. 


10.  SUMMARY 


The  foregoing  review  of  recent  developments  in  radiation-belt 
physics  is  intended  as  a general  overview  of  the  topics  selected.  The 
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literature  contains  many  papers  not  mentioned  here,  and  the  selection 
of  citations  has  admittedly  been  somewhat  parochial.  The  ;L*eader  is 
asked  to  be  tolerant  of  these  shortcomings,  as  the  subject  of  gecmagneti- 
cally  trapped  radiation  is  too  large  in  scope  to  be  covered  thoroughly  in 
the  space  available  here.  However,  it  seems  appropriate  in  closing  to 
mention  once  again  the  outstanding  questions  that  remain  open  in  the 
general  areas  of  investigation  reviewed  above. 

On  the  subject  of  adiabatic  drift  shells,  there  is  still  the  question 
of  whether  quasi -trapping  affects  any  particles  liaving  K = 0 (neglecting 
electrostatic -field  effects)  in  a realistic  model  of  the  outer  magnetosphere, 
i.  e.  , in  a model  whose  surrounding  electrical  currents  close  in  a reason- 
able way.  On  the  topic  of  pitch-angle  diffusion,  one  can  hope  to  include 
self-consistently  the  exchange  of  energy  between  waves  and  particles, 
so  as  to  identify  the  source(s)  of  waves  responsible  for  parasitic  diffusion. 
This  procedure,  and  the  inclusion  of  propagation  effects,  would  enable 
one  to  obtain  the  spectrum,  angular  distribution,  and  spatial  distribution 
of  wave  energy  self-consistently,  without  recourse  to  the  fiat  of  a given 
wave  distribution.  Moreover,  the  inclusion  of  radial  diffusion  and  pitch - 
angle  diffusion  on  an  equal  footing  would  make  it  possible  to  calculate 
both  the  radial  distribution  and  the  pitch-angle  distribution  of  inner- zone 
and  outer-zone  electrons. 

Temporal  variations  of  the  transport  coefficients  represent  the 
major  uncertainty  in  empirical  and  theoretical  studies  involving  radial 
diffusion.  In  the  case  of  outer-zone  protons  and  electrons,  the  problem 
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is  to  account  for  temporal  variations  of  tKe  particle  intensities  in  terms 
of  temporal  variations  of  with  K^.  In  the  case  of  inner -zone 

protons,  the  problem  is  to  model  the  variation  of  transport  coefficients 
and  boundary  conditions  (a)  over  the  solar  cycle,  and  (b)  over  the  quasi- 
periodic  (<*-  8000-yr)  variation  of  the  earth's  magnetic -dipole  moment. 
Temporal  variations  also  play  a complicating  role  in  the  interpretation 
of  data  on  the  alpha/proton  intensity  ratio.  In  this  problem  it  is  certainly 
best  to  view  the  alpha -particle  and  proton  distributions  as  two  separate 
dynamical  entities  occupying  the  same  environment.  Temporal  variations 
of  their  intensity  ratio  thus  emerge  as  incidental  consequences  of  the 
more  interesting  temporal  variations  that  affect  each  particle  species 
separately.  Finally,  to  predict  the  consequences  of  cold-plasma  injec- 
tion into  the  magnetosphere  (an  experimental  project),  one  requires  a 
deeper  understanding  of  the  various  processes  that  control  the  natural 
magnetospheric  environment  and  its  particle  populations. 

These  comments  conclude  the  present  review  of  recent  progress 
in  the  field  of  radiation-belt  phenomenology,  with  emphasis  on  the  past 
four  years  of  research.  Significant  steps  have  been  taken  to  consolidate 
and  refine  a variety  of  ideas  that  had  previously  occupied  the  realm  of 
qualitative  speculation.  Despite  this  progress,  many  of  the  quantitative 
questions  remain  open.  Thus,  in  view  of  the  work  that  must  yet  be  done, 
one  can  only  hope  that  radiation-belt  physics  will  survive  as  an  active 
discipline  for  at  least  "four  more  years". 
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LABORATORY  OPERATIONS 

Th«  Laboratory  Oparatlonr  of  The  Aorotpace  Corporation  ia  conducting 
expai'imantal  and  thaoretlcal  invcstigationa  nacaaaary  for  the  evaluation  and 
application  of  icientific  advancea  to  new  military  concepta  and  ayatema.  Ver- 
aatility  and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
peraomiel  in  dealing  aith  the  many  problema  encountered  in  the  nation'a  rapidly 
developing  apace  and  miaaile  ayatema,  Expertiae  in  the  lateat  aclentific  devel- 
opmenta  ia  vital  to  the  accompliahment  of  taaka  related  to  theae  problema.  The 
laboratorlea  that  contribute  to  thia  reaearch  are: 

Aerophyaica  Laboratory:  Launch  and  reentry  aarodynamica,  haat  trana- 
far,  reentiV  phyalca,  chemical  ki-ietica,  atructural  mechanica,  flight  dynamica, 
atmoapharic  pollution,  and  high-power  gaa  laaera. 

Chemlatry  and  Phyalca  Laboratory:  Atmoapheric  reactiona  and  atmoa- 
pharic  optica,  chemical  reactlona  in  polluted  atmoapherea,  chemical  reactiona 
of  excited  apeclea  in  rocket  plumea,  chemical  thermodynamica.  plaama  and 
laaar-induced  reactiona,  laaer  chemlatry,  propulaion  chemiatry,  apaca  vacuum 
and  radiation  affacta  on  materiala,  lubrication  and  aurface  phenomena,  photo- 
aenaltlva  materiala  And  acnaora,  high  preciaion  laaer  ranging,  and  the  appli- 
cation of  phyalca  and  chemiatry  to  problema  of  law  enforcement  and  biomedicine. 

Electronica  Reaearch  Laboratory;  Electromagnetic  theory,  devicaa,  and 
propagation  phenomena,  including  plaama  electromagnetica;  quantum  electronica, 
laaera,  and  electro-optica;  communication  aciencea,  applied  electronica,  aemi- 
conducting,  auperconducting,  and  cryatal  device  phyaica,  optical  and  acouatical 
imaging;  atmoapheric  pollution;  millii.ieter  wave  and  far-infrared  technology, 

Materiala  Sciencea  Laboratory;  Development  of  new  materiala;  metal 
matrix  compoaitea  and  new  forma  of  carbon;  teat  and  evaluation  of  graphite 
and  ceramic!  in  reentry;  apacecraft  materiala  and  electronic  component!  in 
nuclear  weapon!  environment;  application  of  fracture  mechanica  to  atreaa  cor- 
roaion  and  fatigue-induced  fracture!  in  atructural  metala. 

Space  Phyaica  Laboratory:  Atmoapheric  and  ionoapheric  phyaica,  radia- 
tion from  the  atmoaphere,  denaity  and  compoaition  of  the  atmoaphere,  aurorae 
and  airglow;  magnetoapheric  phyaica,  coamic  raya,  generation  and  propagation 
of  plaama  wavea  in  the  magnetoaphere;  aolar  phyaica,  atudiea  of  aolar  magnetic 
ficlda;  apace  aatronomy,  x-ray  aatronomy;  the  effecta  of  nuclear  exploaiona, 
magnetic  atorma,  and  aolar  activity  on  the  earth'a  atmoaphere,  ionoaphere,  and 
magnetoaphere;  the  effecta  of  optical,  electromagnetic,  and  particulate  radia- 
tiona  in  apace  on  apace  ayatema. 
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